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ABSTRACT 


An  algorithm  and  associated  FORTRAN  77  subroutine,  called  VENTCF2,  is  presented  for  calculating  the 
effects  on  two-layer  compartment  fire  environments  of  the  quasi-steady  flow  through  a circular,  shallow 
(i.e.,  small  ratio  of  depth  to  diameter),  horizontal  vent  connecting  two  spaces.  The  two  spaces  can  be 
either  two  inside  rooms  of  a multi-room  facility  or  one  inside  room  and  the  outside  ambient  environment 
local  to  the  vent.  The  description  of  the  flow  through  the  vent  is  determined  by  combining  considerations 
of  the  uni-directional-type  of  flow  driven  by  a cross-vent  pressure  difference  and,  when  appropriate,  the 
combined  pressure-  and  buoyancy-driven  flows  which  occur  when  the  density  configuration  across  the 
vent  is  unstable,  i.e.,  a relatively  cool,  dense  gas  in  the  upper  space  overlays  a less  dense  gas  in  the 
lower  space.  In  the  algorithm,  calculation  of  the  rates  of  flow  exchange  between  the  two  spaces  is  based 
on  previously  reported  model  equations.  Characteristics  of  the  geometry  and  the  instantaneous 
environments  of  the  two  spaces  are  assumed  to  be  known  and  specified  as  inputs.  The  outputs 
calculated  by  the  algorithm/subroutine  are  the  rates  and  the  properties  of  the  vent  flow  at  the  elevation 
of  the  vent  as  it  enters  the  top  space  from  the  bottom  space  and/or  as  it  enters  the  bottom  space  from 
the  top  space.  Rates  of  mass,  enthalpy,  and  products  of  combustion  extracted  by  the  vent  flows  from 
upper  and  lower  layers  of  inside  room  environments  and  from  outside  ambient  spaces  are  determined 
explicitly.  VENTCF2  is  an  advanced  version  of  the  algorithm/subroutine  VENTCF  in  that  it  includes  an 
improved  theoretical  and  experimental  basis.  The  subroutine  is  completely  modular  and  it  is  suitable  for 
general  use  in  two-layer,  multi-room,  zone-type  fire  model  computer  codes.  It  has  been  tested  over  a wide 
range  of  input  variables  and  these  tests  are  described. 


Keywords: 


building  fires;  compartment  fires;  computer  models;  fire  models;  fire  research; 
mathematical  models;  vents;  zone  models 


INTRODUCTION 


This  work  presents  an  algorithm  and  associated  FORTRAN  77  subroutine,  called  VENTCF2,  to  calculate 
under  arbitrary  conditions  the  instantaneous  effects  on  two-layer  fire  environments  of  the  quasi-steady  flow 
through  a horizontal  vent  connecting  two  spaces  involved  in  a compartment  fire.  The  subroutine  is 
designed  to  be  modular  and  easy  to  integrate  into  any  two-layer  zone-type  compartment  fire  model.  The 
algorithm  VENTCF2,  which  was  developed  in  [2],  is  based  on  ideas  outlined  in  [1],  theoretical  consider- 
ations of  [3],  and  experimental  data  from  reduced-scale  hot-air/cool-air  experiments  of  [4]  and  salt-wa- 
ter/fresh-water experiments  of  [5]  and  [6].  VENTCF2  is  an  advanced  version  of  the  algorithm/subroutine 
VENTCF  [1,  7]  which  was  developed  without  benefit  of  [3]-[5]. 

Flow  through  horizontal  vents  is  a general  problem  associated  with  ventilation  of  enclosed,  heated/cooled 
spaces.  It  is  a problem  whose  general  solution  is  required,  for  example,  if  one  is  to  be  able  to  predict  the 
spread  of  smoke  (i.e.,  fire-heated  and  -contaminated  air)  and  the  flow  of  fresh  air  (i.e.,  oxygen,  which 
could  sustain  a fire,  lack  of  which  could  extinguish  a fire)  during  fires  In  multi-room  facilities.  Reference 
here  is  to  smoke  spread  between  contiguous  rooms,  or  between  a smoky  room  and  the  outside 
environment,  separated  by  a horizontal  partition  (i.e.,  ceiling/floor)  with  penetrations  (i.e.,  vents),  where 
room-to-room  or  room-to-outside,  cross-vent,  pressure  differences  of  arbitrary  magnitude  and  direction 
can  be  generated  by  forced-ventilation  HVAC  systems,  thermal  buoyancy  forces  (i.e.,  stack  effect),  and/or 
wind  effects.  The  problem  has  application  in  fire  scenarios  Involving  top-vented  atria,  stairwells,  ship 
holds,  etc.  The  purpose  for  this  work  Is  to  provide  a computational  tool  that  can  be  used  to  predict  these 
phenomena. 


THE  BASIC  PROBLEM 

Consider  the  flow  through  a horizontal  ceiling/floor  vent  connecting  two  spaces,  one  on  top  of  the  other, 
involved  In  a fire-generated  environment.  The  two  spaces  can  be  two  inside  rooms  of  a multi-room  facility, 
as  depicted  in  Figure  1a  of  APPENDIX  A (page  VENTCF2  - 15),  or  they  can  be  comprised  of  one  inside 
room  of  a facility  and  one  outside  space,  either  above  or  below  the  room,  which  is  used  to  simulate 
conditions,  local  to  the  vent,  of  the  outside  ambient.  The  latter  configurations  are  depicted  in  Figures  1b 
and  1c  of  APPENDIX  A,  respectively. 

When  the  upper  gas  is  less  dense  than  the  lower  gas,  i.e.,  the  fluid  configuration  is  stable,  the  flow 
through  the  vent  is  determined  by  a traditional  orifice-type  flow  model.  Then,  flow  is  determined  by  the 
cross-vent  pressure  difference  without  any  regard  for  buoyancy  effects  (see,  e.g..  References  [8]  and  [9]). 
When  the  configuration  is  unstable  and  the  upper  gas  is  more  dense  than  the  lower  gas,  the  effects  of 
combined  pressure  and  buoyancy  forces  can  be  significant.  For  example  if  the  cross-vent  pressure  is 
relatively  small,  less  than  the  critical  value,  Appi^^oo*  unstable  density  configuration  leads  to  an 
exchange-type  of  flow,  with  gas  in  the  lower  space  rising  into  the  upper  space  and  gas  from  the  upper 
space  dropping  into  the  lower  space,  where  the  flow  rate  from  the  high-  to  the  low-pressure  side  of  the 
vent  is  the  larger  of  the  two.  Even  when  the  cross-vent  pressure  difference  is  large  enough  to  produce 
uni-directional  flow,  the  effect  of  buoyancy  can  be  great  enough  to  reduce  significantly  the  flow  rate  from 
what  it  would  be  in  the  absence  of  a cross-vent  density  difference.  In  the  present  algorithm,  for  the  case 
of  unstable  configurations  the  calculation  of  the  flow  between  the  two  spaces  is  based  on  the  analysis 
of  [2],  which  uses  theoretical  considerations  of  [3],  and  experimental  data  of  [4]  (hot-air/cool-air  in  the  uni- 
directional flow  regime,  flow  from  top  to  bottom)  and  of  [5]  and  [6]  (salt-waterAresh/water  vent  flows  in 
the  exchange-flow  regime). 

The  algorithm/subroutine  VENTCF2  would  be  used  to  calculate  the  net  instantaneous  rates  of  addition  of 
mass,  enthalpy,  and  products  of  combustion  of  interest  and  the  properties  of  the  vent  flows  to  each  of 
the  two  connected  spaces  at  the  elevation  of  the  vent.  A determination  of  where  such  flows  go  once  they 
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enter  the  receiving  spaces  would  be  determined  with  the  use  of  additional  algorithms  and  associated 
subroutines.  These  additional  algorithms,  which  would  use  the  output  of  the  present  algorithm/subroutine, 
would  be  based  on  considerations  beyond  the  scope  of  the  present  work. 

A flow  through  the  horizontal  vent  which  enters  one  of  the  spaces  joined  by  the  vent  is  extracted  from  the 
other  space.  Depending  on  the  configuration  of  the  two  spaces,  the  direction  of  the  flow,  and,  in  the  case 
of  inside  rooms,  the  elevation  of  the  two-layer  interface  (i.e.,  at  the  floor,  ceiling,  or  in-between),  the 
present  algorithm  determines  explicitly  the  rates  of  extraction  of  mass,  enthalpy,  and  products  of 
combustion  from  the  upper  and  lower  layers  of  the  one  or  two  inside  rooms  and/or  from  an  outside 
ambient  space  joined  by  the  vent  under  consideration.  Along  with  other  components  of  flow  to  or  from 
the  layers  of  the  inside  rooms,  determined  with  the  use  of  other  algorithms,  these  rates  would  be  used 
to  continue  in  time  the  solution  to  the  equations  of  the  overall  fire  model.  These  are  the  equations  used 
to  simulate  mathematically  the  facility’s  overall  dynamic  fire  environment. 


LIMITATIONS  ON  VENT  SHAPE 

For  unstable  density  configurations,  the  model  of  [2]  and,  therefore,  the  algorithm/subroutine  presented 
here  is  for  flow  through  a circular,  shallow  (i.e.,  small  ratio  of  depth-to-diameter),  horizontal  vent.  It  is 
expected  that  the  model  will  give  reasonable  estimates  of  flow  even  for  non-circular  vents,  provided  the 
aspect  ratio  (maximum-to-minimum  span)  of  a vent  shape  of  interest  is  not  too  much  different  than  unity. 
Indeed,  one  of  the  example  calculations  of  [2],  which  includes  comparisons  with  some  relevant 
experimental  data,  provides  limited  support  for  the  applicability  of  the  model  in  the  case  of  square  vents. 
However,  when  cross-vent  pressure  differences  are  small-to-moderate  compared  to  AppLooo'  of  the 
model  for  high-aspect-ratio  and/or  moderate-to-large-depth  vents  is  not  valid.  Results  beyond  those 
developed  in  [2]  are  required  before  the  present  work  can  be  evaluated  for  Its  use  In  predicting  flows 
through  the  latter  types  of  vent  shapes. 


THE  ALGORITHMS  AND  ASSOCIATED  FORTRAN  77  SUBROUTINES  VENTCF2  AND  VENTCF2A  FOR 
CALCULATING  THE  EFFECTS  OF  FLOW  THROUGH  HORIZONTAL  CEILING/FLOOR  VENTS 

The  algorithm  VENTCF2  and  a listing  of  its  associated  subroutine,  coded  in  FORTRAN  77,  is  presented 
In  APPENDIX  A.  This  is  a stand-alone  document  that  can  be  Inserted  as  a new  entry  into  the  catalog  of 
modular  algorithms  and  associated  FORTRAN  77  subroutines  of  [10].  APPENDIX  B includes  stand-alone 
documentation  and  the  associated  subroutine  for  VENTCF2A,  a modification  of  VENTCF2.  VENTCF2A 
provides  special  considerations  for  'smoothing'  rates  of  layer  extraction  from  the  flow  source  room  at 
times  of  relatively  thin  adjacent-vent  layers.  When  used  in  a full  zone  model,  and  depending  on  the 
integration  software  for  the  particular  model,  the  considerations  in  VENTCF2A  eliminate  singularities  that 
may  cause  convergence  problems  In  fire  simulations  at  times  when  adjacent-vent  layers  are  growing  or 
shrinking  from  near-zero  depths. 

Reference  [10]  Is  a catalog  of  algorithms/subroutines  useful  for  simulating  the  physical  phenomena  In 
multi-room  zone-type  compartment  fire  model  computer  codes.  Pagination  of  the  two  documents  of  the 
present  APPENDIXES  A and  B,  and  of  all  entries  of  the  Reference  [10]  catalog  is  according  to  name  of 
the  particular  algorithm/subroutine  (in  this  case,  VENTCF2  and  VENTCF2A)  and  page  number  of  the 
catalog  entry. 

The  catalog  of  [10]  was  conceived  of  as  a growing  document  where  the  entries  would  be  available  for 
general  use  by  people  Interested  In:  developing  or  improving,  for  their  own  particular  needs,  a general 
or  special-purpose  multi-room  zone-type  compartment  fire  model;  or  predicting  isolated  compartment  fire 
phenomena,  for  whatever  reason. 
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The  development,  technology  transfer,  and  use  of  a Reference- [10] -type  of  catalog  of  algorithms  and 
associated  subroutines  is  enhanced  by  maintaining  guidelines  for  a uniform  format  of  algorithm/subroutine 
documentation.  In  this  regard  a prototype  format  was  developed  and  used  in  ail  Reference- [10] 
algorithm/subroutine  catalog  entries.  This  format,  which  is  followed  here.  Includes  the  following  elements: 


TITLE  - Should  indicate  the  main  purpose  of  the  algorithm/subroutine. 

DESCRIPTION  - General  description  of  the  algorithm. 


OUTPUT  - List  of  output  variables,  including  definitions  and  units. 

INPUT  - List  of  input  variables,  including  definitions  and  units. 

CALCULATIONS  - Concise  description  of  the  rules  for  obtaining  the  output  variables 

from  the  Input  variables.  This  would  Include  or  refer  explicitly  to 
all  equations  required  in  the  calculation.  If  other  algo- 
rithms/subroutines are  required,  then  these  should  be  readily 
available  and  referenced. 


SUBROUTINES 

USED 


REFERENCES 

SUBROUTINE 

VARIABLES 


PREPARED  BY 


SUBROUTINE 


TESTING  THE  SUBROUTINES 
Subroutine  VENTCF2 

Extensive  parametric  testing  of  the  subroutine  VENTCF2  has  been  carried  out.  A wide  range  of 
environment  scenarios  were  considered  for  each  of  the  three  basic  configurations  of  Figure  1 of  VENTCF2 
(page  VENTCF2  - 15  of  the  APPENDIX  A).  For  each  configuration,  parameters  were  varied  in  a manner 
as  to  simulate  all  possible  combinations  of  the  following: 

a for  the  two  layers  of  an  inside  space:  the  usual  case  with  two  non-zero-thickness  layers 
(i.e.,  layer  interface  between  the  ceiling  and  the  floor),  or  one  non-zero-thickness  layer 
and  one  zero-thickness  layer  (i.e.,  layer  interface  at  the  ceiling  or  floor); 

b.  stable  or  unstable  cross-vent  density  configuration; 


Listing  of  or  explicit  reference  to  each 
algorithm/subroutine  used  to  carry  out  the 
calculations. 

A list  of  references. 

cross-reference  of  ail  nomenclature  (including 

units)  introduced  in  the  above  sections  to  the  nomenclature  used 

In  the  FORTRAN  77  subroutine. 

Names  of  those  who  prepared  the  aigorithm/subroutine  and  date 
of  preparation. 

Listing  of  the  subroutine.  This  would  be  well-commented  and 
would  include  a summary  of  the  purpose  of  the  subroutine  and 
definitions  (including  units)  of  its  input  and  output  variables. 
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c.  the  reference  elevation  for  an  outside  space  is  above,  at,  or  below  the  vent  elevation; 

d.  1.0  atmosphere  reference  pressure,  Pref.v  space  (i.e.,  p^^j  = 101325Pa  with 

<^Prefi  = Pref.1  ■ Pdat  = reference  pressure  in  the  bottom  space,  Pref,2«  varying 

from  0.01  atmospheres  to  2.0  atmospheres  (i.e.,  -O.OOpQ^y  < 5pR£p2  - Pdat) I 

e.  1.0  atmosphere  reference  pressure,  Pref.2'  bottom  space  (i.e.,  p^^^j  = 101325Pa 
with  6pR£p2  = Pref.2  ■ Pdat  = 0)  reference  pressure  in  the  top  space,  Preri,  varying 
from  0.01  atmospheres  to  2.0  atmospheres  (i.e.,  -O.OOPd^t  < <5Pref.i  ^ Pdat) 

Five  distinct  elevations  (relative  to  a datum  elevation),  y^,  N = 1 to  5,  and  three  distinct  densities,  N 
= 1 to  3,  were  required  to  construct  the  environment  scenarios.  These  were  chosen  to  be: 


y^  = 0.0m;  ^2  = Va  = y4  = ^ 

= I.OOkg/m^;  P2  = 0.50kg/m^;  p^  = 0.25kg/m^ 


(1) 


When  the  environment  of  an  inside  room  involved  two  layers,  the  layer  densities  were  specified  to  be  in 
a stable  configuration  (i.e.,  a low  density  layer  above  a high  density  layer). 

All  scenarios  involve  an  example  vent  of  area  \ = 1 m^. 

Table  1 identifies  all  configurations  and  environment  specifications  that  were  used  in  the  parametric 
testing  of  the  subroutine.  This  involved  1 62  basic  Cases.  Each  Case  involved  a separate  parametric 
study  of  over  two  thousand  calls  to  the  subroutine.  These  covered  the  above  pressure  range  of  either 
condition  d.  or  condition  e.  The  calls  involved  incremental  changes  of  the  relative  reference  pressures, 
6Prepi  for  condition  e and  5pp£p2  for  condition  d,  which  were  small  enough  to  reveal  details  of  the 
exchange-flow  phenomena,  these  being  sensitive  to  fractional-pascal-level  variations  in  cross-vent 
pressure  differences. 

For  each  Case,  Table  1 identifies  values  of  the  parameters  which  define  the  configuration  of  the  two 
spaces,  as  illustrated  in  Figure  1 of  VENTCF2,  and  the  states  of  the  environments  in  the  spaces.  These 
are  the  parameters  used  to  call  the  subroutine.  (The  reader  is  referred  to  the  NOMENCLATURE  Section 
for  explanation  of  terms  used  in  the  table.)  Also  explicitly  Identified  In  the  table  are  those  Cases  which 
involve  unstable  cross-vent  density  configurations.  When  cross-vent  pressure  differences  are  small 
enough  in  magnitude,  such  Cases  lead  to  the  exchange-flow  phenomena  which  are  of  particular  interest. 

Presented  In  Figures  1 and  2 are  plots  of  the  calculated  mass  flow  rates  to  upper  room  1,  and 

to  lower  room  2,  MvEhrT.2'  functions  of  the  computed  cross-vent  pressure  difference,  Ap,  for  Cases  28, 
29,  and  30.  As  in  Eq.  (2)  of  VENTCF2,  Ap,  an  output  variable  of  the  subroutine,  is  defined  as 


Ap  = P2  - Pi 


(2) 


where  p^  and  Pg  are  the  pressures  in  the  top  and  bottom  spaces  at  the  elevation  of  the  vent  (see  Figure 
2 of  VENTCF2). 

As  indicated  in  Table  1 , Cases  28,  29,  and  30  involve  two  inside  rooms  where,  as  depicted  in  Figure  1 a 
of  VENTCF2,  both  of  these  have  an  upper  and  a lower  layer.  Case  28  Involves  a neutrally-stable  cross- 
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vent  density  configuration  since  the  density  immediately  above  the  vent  (in  the  lower  layer  of  room  1), 

= /?2  = 0.50kg/m  , Is  identical  to  the  density  Immediately  below  the  vent  (in  the  upper  layer  of  room  2). 
Cases  29  and  30  Involve  an  unstable  cross-vent  density  configuration  since  the  density  immediately  above 
the  vent,  = Pi  = 1.0kg/m^,  Is  greater  than  the  density  immediately  below  the  vent,  p^^  = P2  - 
0.50kg/m^ 

The  scale  of  the  Ap  abscissa  of  Figure  1 is  relatively  coarse,  and  the  plots  reveal  the  calculated  mass  flow 
rates  through  the  vent  as  a result  of  large  cross-vent  pressure  differences,  i.e.,  when  effects  of 
compressibility  become  significant  as  could  be  the  case,  for  example,  for  fires  in  nearly-hermetically-sealed 
facilities.  In  the  figure,  the  sharp  breaks  In  the  Mvent.i  f^vENT.2  Cases  28  and  30  and  for  Case 

29,  respectively,  correspond  to  cross-vent  pressure  differences,  Ap,  which  separate  choked  from 
unchoked  flow  through  the  vent.  Define  Phiqh  ^^icI  Plow  the  absolute  pressures  on  the  high-  and  low- 
pressure  side  of  the  vent,  respectively.  Then,  consistent  with  Eqs.  (7)-(15)  of  VENTCF2,  such  choking 
occurs  for  air  (i.e.,  7 = 1.40)  when  [9,  11] 


Plow/Phigh  S [2/(7  + 1)]^'^  = 0.528 


(3) 


As  indicated  in  Figure  1 , for  Case  29  the  parameters  specified  in  Table  1 lead  to  a prediction  of  choked 
vent  flow  from  room  1 to  room  2 approximately  when  Ap  < -0.472Pd^j  = -0.478(10®) Pa.  Similarly,  for 
Cases  28  and  30  the  parameters  specified  In  Table  1 lead  to  a prediction  of  choked  vent  flow  from  room 
1 to  room  2 approximately  when  Ap  < -0.894p[3^y  = -0.906(1 0®)Pa.  Also,  for  Cases  28  and  30,  choked 
flow  from  room  2 to  room  1 Is  predicted  when  Ap  > = 0.478(10®)  Pa.  Finally,  for  Case  29, 

choked  flow  from  room  2 to  room  1 is  predicted  when  Ap  > 0.894po^y  = 0.906(1 0®)Pa. 

The  scale  of  the  Ap  abscissa  of  Figure  2 is  relatively  fine,  and  the  plots  reveal  the  calculated  mass  flow 
rates  through  the  vent  when  the  cross-vent  pressure  differences  are  at  or  close  to  zero.  Since  the  Case- 
28  density  configuration  is  neutrally-stable,  i.e.,  only  uni-directional  or  zero  flow  is  possible,  the  mass  flow 
rates  through  the  vent  corresponding  to  Ap  = 0 are  seen  to  be  identically  zero.  However,  for  the  unstable 
cross-vent  density  configurations  of  Cases  29  and  30,  the  vent  flow  algorithm  is  seen  to  lead  to  the 
exchange-flow  phenomenon.  In  particular,  there  is  a non-zero  vent  flow  exchange  between  the  top  and 
bottom  spaces  whenever  |Ap|  < AppLooc  where  AppLooD  defined  in  Eq.  (20)  of  VENTCF2. 


Subroutine  VENTCF2A 

The  subroutine  VENTCF2A  was  found  to  give  plausible  results  when  tested  in  CCFM. VENTS  [10, 12-14] 
with  a variety  of  time-dependent  one-  and  two-room  fire  scenarios.  Except  for  the  small-scale,  limited- 
regime  data  of  [4-6],  which  were  used  to  develop  VENTCF2,  the  algorithm  has  not  been  validated 
experimentally. 
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NOMENCLATURE 


Av 

Area  of  vent 

^VENT.I'  1 = 1 to  2 

Mass  flow  rate  of  vent  flow  component  entering  space  I 

Pdat 

Datum  absolute  pressure 

Phigh 

Maximum  of  (p^,  Pg) 

Plow 

Minimum  of  (p^,  Pg) 

Pref.i*  • = Ii  2 

Absolute  pressure  in  space  I at  reference  elevation  yp^p, 

Pi  [Pd 

Absolute  hydrostatic  pressure  in  top  [bottom]  space  at  vent  elevation 

Vc*  Vf  Yr  space  1 

Yceil.i»  Ylayer,!'  Yref.i*  respectively 

VCEIUI  tYREF.ll*  1 = ■•.  2 

If  space  I is  inside  room;  elevation  of  ceiling  [floor]  of  room  I above  datum 
elevation;  if  space  I is  “outside  room':  ycEim  snd  yppp,  are  identical  and  equal  to 
reference  elevation  of  space  I above  datum  elevatiori,  i.e.,  ycEiui  = Yref.i 

VLAYEai’  1 = "It  2 

If  space  I is  inside  room;  elevation  of  upper/lower  layer  interface  in  room  I above 
datum  elevation;  if  space  I is  “outside  room':  yLAVER.!  - Yref.i 

Yvent 

Elevation  of  vent  above  datum  elevation 

Ap 

P2-Pi 

^Pflood 

Minimum  value  of  |Ap|,  in  cases  of  unstable  cross-vent  density  configurations, 
leading  to  uni-directional  vent  flow;  Eq.  (20)  of  APPENDIX  A 

^Pflood.I'  I = I 2 

AppLooD  onset  of  uni-directional  flow  into  space  1 

^pREF.I*  1 = 2 

Pressure  at  reference  elevation,  yppp,,  in  space  1 above  datum  absolute  pressure, 
PoAT>  ^ sp^oe  i is  inside  room,  then  6pppp|  and  yppp,  correspond  to  pressure  and 
elevation,  respectively,  within  the  room  and  at  floor 

P\j  I = 1i  2 

If  space  1 is  inside  room;  density  of  lower  [upper]  layer  in  room  I if  volume  of 
lower  [upper]  layer  is  non-zero  (if  lower  [upper]  layer  volume  is  zero  then 
[Pu.i]  calculation);  if  space  1 is  “outside  room':  is  uniform  density 

there  and  py , s pj^, 
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TABLE  1: 


Configurations/Environments  Used  to  Test  the  VENTCF2  Subroutine 


= 1.00kg/m^  >92  - 0.50kg/m^  > 0.25kg/m^ 

y-  = 0.0m  < y,  = 1.5m  < y,  = 3.0m  < y^  = 4.5m  < yc  = 6.0m 
Pdat  = 1.0132i(10®)Pa;  Av  = im^ 


CASES  1-54:  SPACE  1 (TOP  SPACE)  IS  INSIDE  ROOM,  SPACE  2 (BOTTOM  SPACE)  IS  INSIDE  ROOM  - 
FIGURE  la  of  APPENDIX  A 


CASE 

NO. 


SPACE  SPACE 

1 2 


Lower 

Layer 

yR=y3 

yL=y5 

yc^ys 

1 

2 

3 

4 

5 * 

6 * 


7 

8 

9 

10 
11 
12 


13 

14 

15 

16 
17 


Lower 

Layer 

yR=yi 

yL=y3 

yc=y3 


IWo 

Layers 

yR=yi 

yL=y2 

yc=y3 


Upper 

Layer 

yR=yi 

yL=yi 

yc=y3 


/^L,1 

PL,^~P^~P\^ 
Pl,i=/>1=PU2 
PLA~P‘\^P\-:2rP2 
PlA  ~p^  ^Pl,2~P2 


^Li  ~Pi*^P\i;i-Pi*^P\^-P^ 
~Pz  ^P\J^  “p2  ^Pl,2  ~P^ 
Pla  ~P^^P\i;^’^PL;^^P^ 

PL^  ”P2=^U,2^/’U2“/’l 
PlA  ~P\  >^U,2“/^2^^L,2“^1 
Pla  “^1  '^P\J;z~P2^PL,2~P^ 


PLA~P2'^Pu;2~P^ 

PLA~P2'^Pu;2~P^ 

Pl,i=Pi=^U,2 

PLA~P^-P\i;2 

P\^^-P^>P\i^-p2 

PLA-P^>Pu:2~P2 


^Pref.i  “®*’’Pdat<^Pref4J^Pdat 
^Pref,2“®-»*Pdat’*^^Pref,i  ^Pdat 
^Pref.i  “®'»“Pda^<^Pref,2^Pdat 
^Pref,2'^®*»"Pdat^^Pref,i  ^Pdat 
^Pref.i  =®-»‘Pdat<^Pref,2^Pdat 
^Pref,2=®*»*Pdat'*^^Pref,i  ^Pdat 


^Pref.i  =®*»"Pdat<^Pref,2^Pdat 
^Pref,2=®*>"Pdat^^Pref,i  ^Pdat 
^Pref.i  “®*»'Pdat^^Pref,2^Pdat 
*Pref,2=®*»“Pdat<*Pref,i  ^Pdat 
^Pref.i  =°*»*Pdat^^Pref,2^Pdat 
^Pref.2“°*»'Pdat<^Pref,i  ^Pdat 


^Pref.i  “^•»“Pdat^^Pref,2^Pdat 
^Pref,2”®-’"Pdat<^Pref,i  ^Pdat 

^pREF,1  “^••■PdAT^^PrEF,2^PdAT 

^Pref.2~®*'‘Pdat^*Pref,i  ^Pdat 
^Pref.i  *=®-*"Pdat<*Pref.2^Pdat 
^Pref.2=®-»"Pdat'^^Pref,i  ^Pdat 


* Unstable  cross-vent  density  configuration 
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TABLE  1:  Configurations/Environments  Used  to  Test  the  VENTCF2  Subroutine 

(Cont’d) 

- 1.00kg/m^  > P2-  0.50kg/m^  > p^-  0.25kg/m^ 

= 0.0m  < y2  = 1»5m  < ys  = 3.0m  < y4  = 4.5m  < yg  = 6.0m 
Pdat  = 1.01 325(10®) Pa;  Ay  = Im^ 


CASES  1-54:  SPACE  1 (TOP  SPACE)  iS  INSIDE  ROOM,  SPACE  2 (BOTTOM  SPACE)  IS  INSIDE  ROOM  • 
(CONTD)  FIGURE  1 a of  APPENDIX  A 


CASE  SPACE  SPACE 

NO.  1 2 


19 

20 
21 
22 

23 

24 


IWo 

Layer 

yR=y3 

yL=y4 

yc=ys 


* 

* 


* 

* 


Lower 

Layer 

yR=yi 

yL=y3 

yc=y3 


TWO 

Layers 

yR=yi 

yL=y2 

yc=y3 


Upper 

Layer 

yR=yi 

yL=yi 

yc=y3 


PL^-P^-Pl^ 

PL,^-P^^P^^^P2 

Pu^^P^^PL;^~P2 


Pu\  ~Pz  ^P\i:i  ~P2  <^U2  ~P-\ 
PlA  ^P2 ~P2 ^P\^“P\ 
Pla  ~P2~P\i:i!^Pv:2~Py 
PlA  ”^2”/’u.2^/^U2”/^1 
Pla  ^Pyi^Z^Pi^PLATP^ 
Pla  “^1  ^P\i;2~P2^PL:2~P^ 


Pla  ~P2^Pm^~P'\ 

PLA~P2'^Pyia~P^ 

PLA~P^~Py^^ 

PLA~P^~PM^ 

PLA~Py^Pyi;2~P2 


^Pref.i  =®*»"Pdat'^^Pref,2^Pdat 
^Pref.2=®’*'Pdat'*^*Pref,i  ^Pdat 
*Pref,i  =®’»‘Pdat'*^*Pref.2^Pdat 
^Pref,2”®'»"Pdat<*Pref,i  ^Pdat 
^Pref.i  =®*»“Pdat'^^Pref,2^Pdat 
^Pref,2=°‘»“Pdat'*^^Pref,i  ^Pdat 


^pREF.I  ”®*»“PdAT^^PrEF^^PdAT 

*Pref.2=0**“Pdat<^Pref,i  ^Pdat 

^pREF.I  ■*®**"PdAT<^PrEF,2^PdAT 

*Pref.2=0-*“Pdat<^Pref,i  ^Pdat 
^Pref.1  =0-*-Pdat<^Pref.2^Pdat 
^Pref.2”°-»*Pdat<^Pref,i  ^Pdat 


*Pref,i  =®*»"Pdat<^Pref.2^Pdat 
*Pref.2=0-»"Pdat<^Pref,i  ^Pdat 
^Pref.i  =®*»*Pdat'^^Pref^^Pdat 
^Pref,2“®*»"Pdat^^Pref,i  ^Pdat 
^Pref.i  =0--Pdat<^Pref,2^Pdat 
^Pref.2“®*»"Pdat’*^^Pref,i  ^Pdat 


* Unstable  cross-vent  density  configuration 
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TABLE  1: 
(Cont’d) 


CASES  1-54: 
(CONT’D) 

CASE 

NO. 


* 

* 


* 

* 


* 

* 


Configurations/Environments  Used  to  Test  the  VENTCF2  Subroutine 

= 1.00kg/m^  > P2-  0.50kg/m^  > p^-  0.25kg/m^ 

= 0.0m  < y2  = 1-5m  < ys  = 3.0m  < y4  = 4.5m  6.0m 

Pdat  = 1.01325(10®)Pa;  Ay  = Im^ 


SPACE  1 (TOP  SPACE)  IS  INSIDE  ROOM,  SPACE  2 (BOTTOM  SPACE)  IS  INSIDE  ROOM  - 
FIGURE  la  of  APPENDIX  A 


SPACE  SPACE 

1 2 


Upper 

Layer 

yR=y3 

yL=y3 

yc=ys 


Lower 

Layer 

yR=yi 

yL=y3 

yc=y3 


TWO 

Layers 

yR=yi 

yL=y2 

yc=y3 


Upper 

Layer 

yR=yi 

yL=yi 

yc=y3 


p\i^~P2^PL;z~P^ 
P\iA  ~P2“^P\.^P-\ 
^U,1  ~p\  ~Pl,2 
P\i,^~P^~PU2 
^U,1  ~P^  ^Pl,2~P2 
PUA  ~P\  ^Pl,2~P2 


P\i,^  ~Pz  ^P\ifi.  ~P2  ^PlJsrP^ 
PuA  ~Pz'^P\i;2~P2'^Pu2~P^ 
PUA  ~P2~P\S;t.'^PU2~P\ 

P\}A  ~P2~P\^;2'^P\^~P^ 
p\i,^  ~p^  ^P\i^~P2^PL,2~P^ 
P\iA  ~P^  ^PM^~P2^PL,2^P^ 


P\S,^~P2^P\i^~P^ 
P\i^~P2^Pu;2~P^ 
P\i^~P^~P\i;2 
PUA  ~p\  ^P\iZ 
^U,1  ~p\  ^/^U,2“/^2 
P\i,’\~P\^P\i;2~P2 


^Pref.i  “^•»‘Pdat<^Pref,2^Pdat 
^Pref,2*^®*»"Pda^^^Pref.i  ^Pdat 
*Pref,i  =®*»’Pdat^^Pref,2^Pdat 
^Pref.2”®-'*Pdat^*Pref,i  ^Pdat 
^Pref.i  “®*»"Pdat<^Pref,2^Pdat 
^Pref,2“®-»“Pdat<^Pref,i  ^Pdat 


^Pref.i  =®‘»"Pdat<*Pref,2^Pdat 
^Pref,2=®’»"Pdat’*^^Pref,i  ^Pdat 
^Pref.i  “®*»’Pdat^^Pref,2^Pdat 
^Pref,2=®*»"Pdat^^Pref,i  ^Pdat 
^Pref.i  =®*>’Pdat^^Pref,2^Pdat 
^Pref,2=°-'*Pdat'<^Pref,i  ^Pdat 


^Pref.i  =°-»'Pdat<^Pref,2^Pdat 
^Pref^=®-»*Pdat<^Pref,i  ^Pdat 
^Pref.i  “^••■Pdat^^Pref,2^Pdat 
^Pref,2”®**"Pdat<*Pref,i  ^Pdat 
^Pref.1  ”®*»'Pdat^^Pref,2^Pdat 
^Pref.2“®-»'Pdat^^Pref,i  ^Pdat 


* Unstable  cross-vent  density  configuration 
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TABLE  1:  Configurations/Environments  Used  to  Test  the  VENTCF2  Subroutine 

(Cont’d) 

= 1.00kg/m^  > p2~  0.50kg/m^  > Pz-  0-25kg/m^ 

= 0.0m  < y2  = "l-Sm  < ys  = 3.0m  < y4  = 4.5m  < ys  = 6.0m 
Pdat  = 1-01 325(10®) Pa;  Ay  = Im^ 


CASES  55-108:  SPACE  1 (TOP  SPACE)  IS  OUTSIDE  SPACE,  SPACE  2 (BOTTOM  SPACE)  IS  INSIDE  ROOM  - 
FIGURE  1b  of  APPENDIX  A 


CASE 

SPACE 

SPACE 

NO. 

1 

2 

55 


All 

One 

Layer 

yR=y4 

yL=y4 

yc=y4 


* 

* 


* 

* 


72 


* 

* 


Lower 

Layer 


yR=yi 

yL=y3 

yc=y3 


Two 

Layers 

yR=yi 

yL=y2 

yc=y3 


upper 

Layer 

yR=yi 

yL=yi 

yc=y3 


PL^^P2'^PU1~P^ 

Pu(^P2'^P\^Py 

Pu\^p\^p\^ 

P\^\^P\-P\^ 

Pu^~p^^Pu^~p^ 

^L,1”A>1>/>L,2=P2 


PlA  -Pz<P\i;2^P^<PL;2rP^ 
Pu^  ^Pz^P\i:2^P^'^Pui^P^ 
Pla  ~P2~P\};i'^P\^~P\ 

Pla  ~P2-P\};2^PL;2-P^ 

Pla  ~Py  ^Pu;2~P2^PLA^P^ 
Pla  '^P\y:i~Pi^P\^^Py 


Pu^-Pi^P^^i^Py 

PLA~P2^P\i;2~P‘[ 

PLA-P^~P\i;2 

PLA~P^~P\};2 

PLA~P^^P\};2^P2 

/)m=Pi>/)u^=/>2 


^pREF.I  ='®*»“PdAT'^*PrEF,2^PdAT 

^Pref,2=®*'"Pdat<*Pref,i  ^Pdat 

^pREF.I  “®’*"PdA^^^PrEF,2^PdAT 

^Pref,2=°*»’Pdat^^Pref,i  ^Pdat 
^Pref.i  =®*»"Pdat'<^Pref,2^Pdat 
^Pref,2~®*»"Pdat^^Pref,i  ^Pdat 


^Pref.i  =®-»"Pdat^^Pref,2^Pdat 
^Pref,2“'®*»"Pdat^^Pref,i  ^Pdat 
^Pref.i  ~®*»"Pdat^^Pref,2^Pdat 
^Pref,2”®*»"Pdat<^Pref,i  ^Pdat 
^Pref.i  “®*»“Pdat'^^Pref,2^Pdat 
^Pref,2“°*»‘Pdat'^^Pref,i  ^Pdat 


^pREF.I  =®*»"PdAT’<^PrEF,2^PdAT 

^Pref,2=®'»"Pdat<^Pref,i  ^Pdat 
^Pref.i  ”®**“Pdat^^Pref,2^Pdat 
^Pref,2=®*»“Pdat^^Pref.i  ^Pdat 

^pREF.I  =®*»"PdAT^^PrEF,2^PdAT 

^Pref,2=®*»'Pdat'*^^Pref,i  ^Pdat 


* Unstable  cross-vent  density  configuration 
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TABLE  1:  Configurations/Environments  Used  to  Test  the  VENTCF2  Subroutine 

(Cont’d) 

= 1.00kg/m^  > P2-  0.50kg/m^  > P3  = 0.25kg/m^ 

= 0.0m  < y2  = 1-5m  < ya  = 3.0m  < y4  = 4.5m  < ys  = 6.0m 
Pdat  = 1.01 325(10®)  Pa;  Ay  = 1m^ 


CASES  55-108:  SPACE  1 (TOP  SPACE)  IS  OUTSIDE  SPACE,  SPACE  2 (BOTTOM  SPACE)  IS  INSIDE  ROOM  - 
(CONT’D)  FIGURE  1b  of  APPENDIX  A 


CASE  SPACE  SPACE 

NO.  1 2 


73 

74 

75 

76 

77 

78 


All 

One 

Layer 

yR=y3 

yL=y3 

yc=y3 


79 

80 
81 
82 

83 

84 


85 

86 

87 

88 

89 

90 


Lower 

Layer 


yR=yi 

yL=y3 

yc=y3 


Two 

Layers 

yR=yi 

yL=y2 

yc=y3 


Upper 

Layer 

yR=yi 

yL=yi 

yc=y3 


P\^^-P2'^P\^-P\ 
PL,^  ~p^  ~Pl^ 
PL^-P^-P\^2 
/>L,1=A>1>/>L^=P2 
PU^-P^^P\^^~P2 


PL,^  =A>3<Pu,2“^2<^L^”/’i 
Pla  ~Pi^P\i;2~P2’*^P\^~P\ 
Pl,^  ~P2~P\i;z^PL:z~PA 
Pla  -P2~P\s^^PvjrP\ 

PL,^  -P^  ^P^)^^P2^PL^~P^ 
PL,^  ~P-\  ^PU^~P2'^PL^~P^ 


^L,i  ~P2'^Pu^~Py 

PlA  ~P2“^P\i,2~P\ 

PLA-P\-P\i^ 

PL,^-P^-P\i;2 

A>L,i=Pi>^U,2“^2 

PU^~P\^P\i^~P2 


^Pref.i  “®’»"Pdat^^Pref,2^Pdat 
^Pref,2“®*»"Pdat^^Pref,i  ^Pdat 
*Pref,i  =®*»"Pdat<^Pref,2^Pdat 
^Pref,2“®’*’Pdat^^Pref,i  ^Pdat 
^Pref.i  “®*»'Pdat^^Pref,2^Pdat 
^Pref,2  '"Pdat  ^^Pref.i^  Pdat 


^Pref.i  ”®*»"Pdat<^Pref,2^Pdat 
^Pref,2“®**"Pdat^^Pref,i  ^Pdat 
^Pref.i  ”®*»"Pdat^^Pref,2^Pdat 
^Pref,2“®*'"Pdat<^Pref,i  ^Pdat 
^Pref.i  “®*'"PdAT^^PrEF,2^PdAT 
^Pref.2”®‘»"Pdat^^Pref,i  ^Pdat 


^pREF.I  ”°*'"PdAT^^PrEF,2^PdAT 

^Pref,2”®*'*Pdat'*^^Pref,i  ^Pdat 
*Pref,i  “®*'"PdAT^^PrEF,2^PdAT 
^Pref,2”®’'"Pdat^^Pref,i  ^Pdat 
^Pref.1  “®**"Pdat^^Pref,2^Pdat 
^Pref.2”®*'"Pdat^^Pref,i  ^Pdat 


* Unstable  cross-vent  density  configuration 
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TABLE  1:  Configurations/Environments  Used  to  Test  the  VENTCF2  Subroutine 

(Cont’d) 

= 1.00kg/m^  > Pz-  0.50kg/m^  > P3  = 0.25kg/m^ 

= 0.0m  < y2  = l-Sm  < ys  = 3.0m  < y4  = 4.5m  < ys  = 6.0m 
Pdat  = 1.01 325(10®) Pa;  Ay  = Im^ 


CASES  55-108:  SPACE  1 (TOP  SPACE)  IS  OUTSIDE  SPACE,  SPACE  2 (BOTTOM  SPACE)  IS  INSIDE  ROOM  - 
(CONT D)  FIGURE  1 b of  APPENDIX  A 


CASE  SPACE  SPACE 

NO.  1 2 


91 

92 

93 

94 

95 

96 


All 

One 

Layer 

yR=y2 

yL=y2 

yc=y2 


97 

98 

99 

100 
101 
102 


103 

104 

105 

106 

107 

108 


Lower 

Layer 


yR=yi 

yL=y3 

yc=y3 


Two 

Layers 

yR=yi 

yL=y2 

yc=y3 


upper 

Layer 

yR=yi 

yL=yi 

yc=y3 


pU^^P2^PL;2^P^ 

P\i^~P2"^PL,^~P^ 

Py^^-P^-Pv^ 

^U,1”^1“^L,2 

P\},‘\^P\>P\^^P2 

P\i^~P^^PU^'^P2 


PUA  '^P2^P\S;2~P2'^PL^~P^ 
PuA  ~Pz'*^P\^;i~P2’^P\^^P\ 
P\}A  ^P2~P\i;2^Pu2~P^ 

P\i^  -P2^Pm:3l^P\^~P\ 

^U,1  ~P^  '^P\i:t^P2'^PU‘^P^ 
^u,i  ^Py  '^P\}^~P2^PL;2~Py 


P\},y~P2'^P\i;2-Py 
P\iA  ~P2^P\i;2~Py 
P\i,y  ~Py  ~P\i;2 
P\},y~Py~P\i;t 
P\iA~Py^P\i:2~P2 
P\i,\^Py^P\sX^P2 


^PREF.I  ^^’’'PdAT^^PrIF^^PdAT 
*Pref,2=®*»"Pda^'^^Pref,i  ^Pdat 
^Pref.i  ~®'»"Pdat^^Pref.2^Pdat 
*Pref4!=®*»"Pdat^^Pref,i  ^Pdat 
^Pref.i  =®*'"Pdat'^*Pref.2^Pdat 
^Pref.2=®*»"Pdat^*Pref,i  ^Pdat 


^pREF.I  ”°*»"PdAT^^PrEF,2^PdAT 

^Pref^=®*»"Pdat^^Pref,i  ^Pdat 
^Pref.i  ®'®*»"Pdat^^Pref,2^Pdat 
^Pref.2“®*»"Pdat<^Pref,i  ^Pdat 
^Pref,i  =®*»“Pdat'*^^Pref^^Pdat 
^Pref.2=°*»’Pdat'*^^Pref,i  ^Poat 


^Pref.1  =0-»-Pdat<^Pref.2^Pdat 
^Pref.2“®*'"Pdat^^Pref,i  ^Pdat 
^Pref.1  =0--Pdat<*Pref.2^Pdat 
^Pref.2-0-»'Pdat<^Pref,i  ^Pdat 
^Pref.i  “°’»"Pdat<^Pref.2^Pdat 
^Pref.2”®*»"Pdat^^Pref,i  ^Pdat 


* Unstable  cross-vent  density  configuration 
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TABLE  1 : Configurations/Environments  Used  to  Test  the  VENTCF2  Subroutine 

(Cont’d) 

- 1.00kg/m^  > Pz-  0-50kg/m^  > Pz-  0.25kg/m^ 
y.  = 0.0m  < Vo  = 1.5m  < y,  = 3.0m  < y^  = 4.5m  < = 6.0m 

= 1.0132i(105)Pa;  = Im^ 


CASES  109-162:  SPACE  1 (TOP  SPACE)  IS  INSIDE  ROOM,  SPACE  2 (BOTTOM  SPACE)  IS  OUTSIDE 

SPACE  - FIGURE  1c  of  APPENDIX  A 


CASE 

NO. 


SPACE  SPACE 

1 2 


Lower 

Layer 


109 

110 
111 
112 

113  * 

114  * 


yR=y3 

yL=y5 

yc=y5 


All 

One 

Layer 

yR=y4 

yL=y4 

yc=y4 

PL,^~P2'^PL^~P^ 
PL,^~P^~PL,2 
^L,1  “^1  ~Pl,2 
PLA~Py^PL^~P2 
PW\~P\^P\^~P2 


^Pref,i  =®*»’Pdat^*Pref,2^Pdat 
^Pref.2=®*»’Pdat<^Pref,i  ^Pdat 
^Pref.i  ”®*»’Pda^<^Pref4j^Pdat 
^Pref,2~®*»“Pdat<^Pref,i  ^Pdat 
^Pref.i  =®’»'Pdat<^Pref,2^Pdat 
^Pref,2“®‘»“Pdat<^Pref,i  ^Pdat 


115 

116 

117 

118 

119 

120 


yR=y3 

yL=y3 

yc=y3 


^Li  ~Pi^P\};i^P2'^P\^~P\ 
Pl^  “Pz  ^P\i;2  ~P2  ~P^ 

^L,1  ~p2  ”^U,2  ^Puz 
^L1  ~P2~Pvia^PL^~P^ 

Pu^  ~P^  '^P\i;l-P2'^Pu^~P^ 
PL,^  ~p^  '^P\i;2-p2'^PL^~P^ 


^Pref.i  =®*»’Pdat^^Pref,2^Pdat 
*Pref,2=®*»“Pdat<*Pref,i  ^Pdat 
^Pref.i  “®*»’Pdat<^Pref,2^Pdat 
^Pref,2=®**“Pdat^^Pref,i  ^Pdat 
^Pref.i  =®*'’Pdat*^^Pref,2^Pdat 
^Pref,2”®**"Pdat^^Pref,i  ^Pdat 


121 

122 

123 

124 

125  * 

126  * 


yR=y2 

yL=y2 

yc=y2 


^Ui  ^P2'^P\i;irP^ 
Pu^^p2'^P\^:^~P^ 
/>L,1=^1=A>U,2 

^L,1=Pi>A>U.2”^2 

PU‘\-Py^P\i;i-Pz 


^Pref.i  “®**"Pdat^*Pref,2^Pdat 
^Pref,2=®-»"Pdat^^Pref,i  ^Pdat 
^Pref.i  ”®'»"Pdat<^Pref,2^Pdat 
^Pref,2“®**'Pdat'*^^Pref,i  ^Pdat 
^Pref.i  =®-»"Pdat<^Pref.2^Pdat 
^Pref,2”®*»"Pdat^^Pref,i  ^Pdat 


* Unstable  cross-vent  density  configuration 
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TABLE  1:  Configurations/Environments  Used  to  Test  the  VENTCF2  Subroutine 

(Cont’d) 

py  - 1.00kg/m^  > p2~  0-50kg/m^  > Pz-  0.25kg/m^ 

y.  = 0.0m  < y,  = 1.5m  < y,  = 3.0m  < y^  = 4.5m  < w^-  6.0m 

Pdat  = '‘•0132l(10^Pa;  a/=  1m2 


CASES  109-162:  SPACE  1 (TOP  SPACE)  IS  INSIDE  ROOM,  SPACE  2 (BOTTOM  SPACE)  IS  OUTSIDE 

(CONTD)  SPACE  - FIGURE  1c  of  APPENDIX  A 


CASE  SPACE  SPACE 

NO.  1 2 


Two 

Layer 

yR=y3 

yL=y4 

yc=y5 

127 

128 

129 

130 

131  * 

132  * 


All 

One 

Layer 

yR=y4 

yL=y4 

yc=y4 

PL,^  ~P2^PU2~P'y 

Pla  ^Pi^P\^-P\ 

PU-\~P\-P\^ 

PU\-P\-pLa. 

PU\-P\^P\^-P2 

P\^\^P^^P\^-P2 


*Pref,i  ='®’’"Pdat'<^Pref.2“Pdat 
^Pref.2=®-»’Pdat'*^^Pref,i  ^Pdat 
^Pref.i  ”®-»"Pdat^^Pref,2^Pdat 
^Pref,2~®*»"Pdat^^Pref,i  ^Pdat 

^pREF.I  =®’»*PdAT^^PrEF.2^PdAT 

*Pref.2=0-»-Pdat<^Pref,i  ^Pdat 


133 

134 

135 

136 

137  * 

138  * 


yR=y3 

yL=y3 

yc=y3 

Pl^  -Pi^P\i;2~P2'^Pu2~P^ 
PlA  ~Pz^P\};2~P2^PL;2^P^ 
Pla  -P2~P\i;2'^P\^^P^ 

PlA  “^2”^U,2’*^^U2”^1 
PlA  ~P^  ^^U,2“/^2^PU2”/’l 

' ^L,1”^1^^U.2=^2'^^L2“^1 


^PREF,1  =0*»'PdAT<^PrEF,2^PdAT 
^Pref,2=®'»“Pdat^^Pref,i  ^Pdat 
^Pref.i  “®*»‘Pdat'*^^Pref,2^Pdat 
^Pref4S“°*»"Pdat'^^Pref,i  ^Pdat 
^Pref.1  =0-»-Pdat<^Pref,2=sPdat 
^Pref.2=‘®-»*Pdat'<^Pref,i  ^Pdat 


139 

140 

141 

142 

143 

144 


yR=y2 

yL=y2 

yc=y2 

PLA~P2’^P\i;2~^P^ 

PLA~P2"^P\i^^P^ 

PiA~P^~Pyy^ 

PLA^P^~Pu;2 

PLA-P^'^P\J;2-P2 

Pla  ^P\i;2^P2 


^Pref.i  =®*»"Pdat^^Pref.2^Pdat 
^Pref.2“®*»“Pdat^^Pref,i  ^Pdat 
^Pref.i  =®**“Pdat^^Pref,2^Pdat 
^Pref.2=°*»"Pdat^^Pref,i  ^Pdat 
^Pref.1  =0--Pdat<^Pref.2^Pdat 
^Pref,2“°*»’Pdat^*Pref,i  ^Pdat 


* Unstable  cross-vent  density  configuration 
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TABLE  1:  Configurations/Environments  Used  to  Test  the  VENTCF2  Subroutine 

(Cont’d) 

= 1.00kg/m^  > P2  ” 0.50kg/m^  > P3  = 0.25kg/m^ 

= 0.0m  < Y2~  < Va  = y4  ” < Vs  “ 

= 1.01 325(10®) Pa;  = Im^ 


CASES  109-162: 

SPACE  1 (TOP  SPACE)  IS  INSIDE  ROOM,  SPACE  2 (BOTTOM  SPACE)  IS  OUTSIDE 

(CONT’D) 

SPACE  - 

FIGURE  1c  of  APPENDIX  A 

CASE 

SPACE 

SPACE 

NO. 

1 

2 

Upper 

All 

Layer 

One 

Layer 

Yr=Y4 

Yl^Ys 

Yl=Y4 

Yc=Y5 

Yc=Y4 

145 

PV,^~P2^PL,2~P^ 

^PREF.I  “®*'"PdAT<^PrEF,2^PdAT 

146 

Pii,^~P2^PL,2~P^ 

^Pref,2“®*'"Pda^^*Pref,i  ^Pdat 

147 

PU,^~P^~PL^ 

^Pref.1  ”®-*’Pdat<*Pref.2^Pdat 

148 

P^J,^~P^~PL,2 

^Pref,2”®'*’Pdat^^Pref,i  ^Pdat 

149 

* 

P[J,^~P^^PL,2~P2 

^Pref,i  “°-»‘Pdat^^Pref,2^Pdat 

150 

* 

P\iA  ^PL,'2rP2 

^Pref,2”®-*'Pdat^^Pref,i  ^Pdat 

Yr=Y3 

Yl=Y3 

Yc=Y3 

151 

PyiA  -P2^Py^^~P2^Pl^~P^ 

^Pref.i  “°*»*Pdat<^Pref.2^Pdat 

152 

P\iA  ~Pin^P\ip~P9^P\-P“P-\ 

^Pref.2“®*'"Pdat^^Pref,i  ^Pdat 

153 

P\iA  -p2~P\ip'^P\^~P^ 

*pREF,1  “®*»'PdAT<^PrEF,2^PdAT 

154 

^Pref,2“®’'’Pdat<^Pref,i  ^Pdat 

155 

* 

^U,1  ~P\  ^P\^P.~P2'^P\^~P^ 

^Pref.1  =®‘»’Pdat'*^^Pref,2^Pdat 

156 

* 

/^U,1  ~P^  '^P\iprP2'^P\^~P\ 

^Pref,2=®*»*Pdat^^Pref,i  ^Pdat 

Yr=Y2 

Yl=Y2 

Yc=Y2 

157 

P\ip~P2“*^P\iP~P^ 

^Pref.i  =®-»*Pdat^^Pref.2^Pdat 

158 

^Pref.2“®*'"Pdat<^Pref,i  ^Pdat 

159 

^u,i  ^P^  ~Pup. 

^Pref.i  “®**"Pdat<^Pref.2^Pdat 

160 

Pyip“P\~Pyip. 

^Pref.2“®-''Pdat<*Pref,i  ^Pdat 

161 

* 

P\JA  ~P^  '^P\ip-P2 

^Pref.i  “^••"Pdat<^Pref,2^Pdat 

162 

* 

P\i^~P^^P\ip~P2 

^Pref.2"®’»’Pdat^^Pref,i  ^Pdat 

* Unstable  cross-vent  density  configuration 
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VENT, 2 IVIVENT.I 


Ap  [Pa] 


Figure  1 . Plots  of  the  calculated  mass  flow  rate  to  upper  space  1 , and  to  lower  space 

2,  Mvent,2»  a®  functions  of  the  computed  (coarse-scale)  cross-vent  pressure 
difference,  Ap,  for  Cases  28,  29,  and  30  (see  Table  1). 
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Ap  [Pa] 


Figure  2.  Plots  of  the  calculated  mass  flow  rate  to  upper  space  1,  and  to  lower  space 

2,  as  functions  of  the  computed  (fine-scale)  cross-vent  pressure  difference, 

Ap,  for  uases  28,  29,  and  30  (see  Table  1). 
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APPENDIX  A (pages  VENTCF2-1  to  VENTCF2-22) 


VENTCF2  - 


VENTCF2A  - 


CALCULATION  OF  THE  FLOW  THROUGH  A HORIZONTAL  CEILING/FLOOR  VENT 
CONNECTING  TWO  SPACES 


APPENDIX  B (pages  VENTCF2A-1  to  VENTCF2A-25) 


CALCULATION  OF  THE  FLOW  THROUGH  A HORIZONTAL  CEILING/FLOOR  VENT 
CONNECTING  TWO  SPACES  WITH  "SMOOTHING"  OF  LAYER  EXTRACTION  RATES 
AT  TIMES  OF  RELATIVELY  THIN  ADJACENT-VENT  LAYERS 
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VENTCF2  - CALCULATION  OF  THE  FLOW  THROUGH  A HORIZONTAL  CEILING/FLOOR  VENT 
CONNECTING  TWO  SPACES 


DESCRIPTION 

Consider  an  instant  of  time  during  the  simulation  of  a multi-room  compartment  fire  environment.  This 
algorithm  calculates  the  flow  of  mass,  enthalpy,  oxygen,  and  other  products  of  combustion  through  a 
horizontal  vent  located  in  a ceiling/floor  partition  common  to  any  two  inside  rooms  of  the  facility  or 
between  an  inside  room  and  the  outside  environment  local  to  the  vent.  VENTCF2,  is  an  advanced  version 
of  the  algorithm/subroutine  VENTCF  [1]  in  that  it  includes  an  improved  theoretical  [2,  3]  and  experimental 
basis. 

Depicted  in  Figure  1 a is  the  vent  and  the  two  spaces  when  they  are  both  inside  rooms  of  a multi-room 
facility.  Figures  1 b and  1 c depict  the  situation  when  the  two  spaces  involve  one  inside  room  of  the  facility 
and  one  outside  space,  either  above  or  below  the  room,  in  which  Is  simulated  the  outside  environment 
local  to  the  vent. 

As  in  Figure  1 , designate  the  top  space  as  space  1 and  the  bottom  space  as  space  2.  It  is  assumed  that 
the  temperature,  density,  concentration  of  oxygen  and  of  other  products  of  combustion  of  interest  in  the 
upper  and  lower  layer  of  each  inside  room  and  in  the  environment  local  to  the  vent  of  an  outside  space 
are  specified.  Also  specified  In  each  inside  room  are:  the  elevation  above  the  datum  elevation  of  the  floor, 
and  the  upper-layer/lower-layer  interface;  and  the  pressure  at  the  floor  above  the  specified  datum 
pressure.  Specified  in  an  outside  space  are;  a reference  elevation  above  the  datum  elevation,  and  the 
pressure  at  this  reference  elevation  above  the  specified  datum  pressure. 

When  the  upper  gas  is  less  dense  than  the  lower  gas,  i.e.,  the  fluid  configuration  is  stable,  the  flow 
through  the  vent  is  determined  by  a traditional  orifice-type  flow  model.  Then,  flow  is  determined  by  the 
cross-vent  pressure  difference  without  any  regard  for  buoyancy  effects  (see,  e.g.,  [4]  and  [5])  and  the 
present  algorithm/subroutine  is  identical  to  that  of  VENTCF  [1].  When  the  configuration  Is  unstable  and 
the  upper  gas  is  more  dense  than  the  lower  gas,  the  effects  of  combined  pressure  and  buoyancy  forces 
can  be  significant.  For  example  if  the  cross-vent  pressure  is  relatively  small,  the  unstable  density 
configuration  leads  to  an  exchange-type  of  flow,  with  gas  in  the  lower  space  rising  into  the  upper  space 
and  gas  from  the  upper  space  dropping  into  the  lower  space,  where  the  flow  rate  from  the  high-  to  the 
low-pressure  side  of  the  vent  Is  the  larger  of  the  two.  Also,  even  when  the  cross-vent  pressure  difference 
large  enough  to  produce  uni-directional  flow,  the  effect  of  buoyancy  can  be  great  enough  to  reduce 
significantly  the  flow  rate  from  what  it  would  be  in  the  absence  of  a cross-vent  density  difference.  In  the 
present  algorithm,  for  the  case  of  unstable  configurations  the  calculation  of  the  flow  between  the  two 
spaces  is  based  on  the  theory  and  analysis  of  [3]  and  [2],  respectively. 

For  unstable  density  configurations,  the  model  of  [2]  and,  therefore,  the  VENTCF2  algorithm/subroutine 
itself  is  for  flow  through  a circular,  shallow  (i.e.,  small  ratio  of  depth  to  diameter),  horizontal  vent.  It  is  also 
expected  that  the  model  will  give  reasonable  estimates  of  flow  even  for  non-circular  vents,  provided  the 
aspect  ratio  (maximum-to-minimum  span)  of  a vent  shape  of  interest  is  not  too  much  different  than  1 . 
Indeed,  one  of  the  example  calculations  of  [2],  which  includes  comparisons  with  some  relevant 
experimental  data,  provides  limited  support  for  the  applicability  of  the  model  in  the  case  of  square  vents. 
However,  USE  OF  VENTCF2  IN  HIGH-ASPECT-RATIO  AND/OR  MODERATE-TO-LARGE-DEPTH  VENT 
SCENARIOS  IS  NOT  VALID  when  cross-vent  pressure  differences  are  small-to-moderate  compared  to 
AppLooD*  Results  beyond  those  developed  In  [2]  are  required  before  the  present  work  can  be  used  for 
the  latter  types  of  vent  shape. 
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The  geometry  and  the  conditions  local  to  the  vent  which  determine  the  characteristics  of  the  vent  flow  are 
depicted  in  Figure  2.  These  include:  the  densities,  and  P2*  ^nd  the  hydrostatic  pressures,  p^  and  P2, 
at  the  elevation,  but  away  from  the  immediate  vicinity  of  the  vent  in  the  upper  and  lower  spaces, 
respectively,  and  the  area,  A^,  and  shape  of  the  vent.  Regarding  the  shape,  at  the  present  time  results 
for  horizontal  vent  flows  are  only  available  for  circular  or  square  vents.  Other  properties  local  to  the  vent 
and  indicated  in  Figure  2 are  T^  and  T2,  the  absolute  temperatures,  Cq2.i  and  Cq2.2>  concentrations  of 
oxygen,  and  c^^^  and  c^2'  K = 2 to  Mppo^,  the  concentrations  of  a product  of  combustion  K. 

To  simulate  fire  scenarios  at  times  when  vent  flows  are  driven  by  arbitrarily  high  cross-vent  pressure 
differences  (i.e.,  when  compressibility  effects  begin  to  be  significant),  whether  for  stable  or  unstable  cross- 
vent density  configurations,  VENTCF2  implements  the  Ideas  Introduced  in  reference  [5]  and  implemented 
previously  in  the  algorithm/subroutines  VENTHP  [6]  and  VENTCF  [1].  Such  high  cross-vent  pressure 
differences  could  occur,  for  example,  in  fire  scenarios  involving  flows  through  cracks  in  otherwise 
hermetically-sealed  fire  compartments. 

The  VENTCF2  subroutine  has  been  subjected  to  extensive  parametric  testing  In  a "stand-alone"  mode  [7]. 


OUTPUT 

^02,1  1^02.2! 

Concentration  of  oxygen  in  the  top  [bottom]  space  at  the  elevation  of  the  vent,  [(kg  of  oxygen)/(kg 
of  layer)] 


^K.1  I^K.2l 

Concentration  of  product  K in  the  top  [bottom]  space  at  the  elevation  of  the  vent,  [(unit  of  product 
K)/(kg  of  layer)] 

^NT.02.1*  I = 1 or  2 

Concentration  of  oxygen  in  the  component  of  the  vent  flow  entering  space  I,  provided  such  vent 
flow  component  is  non-zero,  i.e.,  provided  is  non-zero,  [(kg  of  oxygen)/(kg  of  vent  flow)] 


^MT.K.p  I = 1 or  2;  K = 2 to  NPROD 


Concentration  of  product  of  combustion  K in  the  component  of  the  vent  flow  entering  space  I, 
provided  such  vent  flow  component  is  nonzero,  i.e.,  provided  is  non-zero,  [(unit  of 

product)/(kg  of  vent  flow)] 

^u.i  I = 1 or  2 

If  space  I is  an  inside  room:  Rate  at  which  mass  is  added  to  the  upper  [lower]  layer  of  room  I 
due  to  the  vent  flow  component  which  enters  the  other  space.  Note  that  this  will  always  be 
negative,  since  the  layer  which  supplies  the  material  flowing  to  the  other  space  will  always  have 
mass  extracted  from  it.  [W] 

If  space  I is  an  "outside  room":  , is  the  rate  at  which  mass  is  added  to  space  I due  to  the  vent 
flow  component  which  enters  the  other  space  (i.e.,  the  inside  room).  Note  that  this  will  always 
be  negative,  since  an  outside  space  which  supplies  material  flowing  through  the  vent  to  an 
adjacent  space  will  always  have  mass  extracted  from  it.  is  identical  to  M^  ,.  [W] 
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^VENT.I’  I - 1 to  2 


Mass  flow  rate  of  the  vent  flow  component  entering  space  I.  [kg/s] 

^02.LI  t^02.u.j  I I = 1 or  2 

If  space  I is  an  inside  room:  Rate  at  which  oxygen  is  added  to  the  upper  [lower]  layer  of  room 
I due  to  the  vent  flow  component  which  enters  the  other  space.  Note  that  this  will  always  be 
negative,  since  the  layer  which  supplies  the  material  flowing  through  the  vent  to  the  other  space 
will  always  have  Its  convected  oxygen  extracted  from  it.  [(kg  of  02)/s] 

If  space  I Is  an  'outside  room':  f^02,u.i  which  oxygen  is  added  to  space  I due  to  the 

vent  flow  component  which  enters  the  other  space  (i.e.,  the  inside  room).  Note  that  this  will 
always  be  negative,  since  an  outside  space  which  supplies  material  flowing  through  the  vent  to 
an  adjacent  space  will  always  have  its  convected  oxygen  extracted  from  it.  I^02.l,i  identical  to 

^ 02,U,I'  02)/S] 

^02,VENT.I'  i = 1 to  2 

Mass  flow  rate  of  oxygen  of  the  vent  flow  component  entering  space  I.  [(kg  of  oxygen)/s] 

^K.ui  I^K.u.|]'  K = 2 to  NppQ^;  I = 1 or  2 

If  space  I is  an  inside  room:  Rate  at  which  product  of  combustion  K is  added  to  the  upper  [lower] 
layer  of  room  I due  to  the  vent  flow  component  which  enters  the  other  space.  Note  that  this  will 
always  be  negative  since  the  layer  which  supplies  the  material  flowing  through  the  vent  to  the 
other  space  will  always  have  Its  convected  product  of  combustion  K extracted  from  it.  [(unit  of 
product  K)/s] 

If  space  I Is  an  'outside  room':  , is  the  rate  at  which  product  of  combustion  K is  added  to 

space  I due  to  the  vent  flow  component  which  enters  the  other  space  (i.e.,  the  inside  room).  Note 
that  this  will  always  be  negative,  since  an  outside  space  which  supplies  material  flowing  through 
the  vent  to  an  adjacent  space  will  always  have  the  convected  product  of  combustion  K extracted 
from  it.  is  identical  to  P,^y [(unit  of  product  K)/s] 

^KVENT.i*  I = 1 to  2;  K = 2 to  Np^^Q 

Flow  rate  of  product  of  combustion  K in  the  vent  flow  component  entering  space  I.  [(unit  of 
product  K)/s] 


Pi  [P2I 


Absolute  hydrostatic  pressure  in  the  top  [bottom]  space  at  the  elevation  of  the  vent.  [Pa  = 
kg/(m  • s^] 

^u.i  I = 1 or  2 

If  space  I is  an  inside  room:  Rate  at  which  enthalpy  added  to  the  upper  [lower]  layer  of  room  I 
due  to  the  vent  flow  component  which  enters  the  other  space.  Note  that  this  will  always  be 
negative,  since  the  layer  which  supplies  the  material  which  flows  to  the  other  space  will  always 
have  its  convected  enthalpy  extracted  from  It.  The  enthalpy  Is  based  on  the  absolute  temperature 
of  the  flow,  Tvent.i-  [W] 
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If  space  I is  an  •outside  room':  Oy , is  the  rate  at  which  enthalpy  is  added  to  space  I due  to  the 
vent  flow  component  which  enters  the  other  space  (i.e.,  the  inside  room).  Note  that  this  will 
always  be  negative,  since  an  outside  space  which  supplies  material  which  flows  through  the  vent 
to  an  adjacent  space  will  always  have  its  convected  enthalpy  extracted  from  It.  The  enthalpy  Is 
based  on  the  absolute  temperature  of  the  flow,  Ty ,.  Is  identical  to  Oy ,.  [W] 


^vENT.i'  I - 1 to  2 


Total  enthalpy  flow  rate  in  the  vent  flow  component  entering  space  I.  This  is  based  on  the 
absolute  temperature  of  the  flow,  Ty^^,.  [W] 

T,F2] 

Absolute  temperature  in  the  top  [bottom]  space  at  the  elevation  of  the  vent.  [K] 

'^'vENT.i’  I = 1 to  2 

Absolute  temperature  of  the  vent  flow  component  entering  space  I,  provided  such  vent  flow 
component  is  non-zero.  [K] 

Ap 

P2  - p^,  i.e.,  difference  between  the  pressure  in  the  bottom  space  at  the  elevation  of  the  vent,  and 
the  pressure  in  the  top  space  at  the  elevation  of  the  vent.  [Pa  = kg/(m  • s^j 

Pi  [P2] 

Density  in  the  top  [bottom]  space  at  the  elevation  of  the  vent,  [kg/mg] 

PvENT.i*  J = "I  to  2 

Density  of  the  vent  flow  component  entering  space  I,  provided  such  vent  flow  component  is 
nonzero,  [kg/m^] 


INPUT 

Ay 

Area  of  the  vent  [m^]. 

^UKi  K = 2 to  NppQQi  I = 1 or  2 

if  space  I is  an  inside  room:  Concentration  of  product  of  combustion  K in  lower  [upper]  layer  of 
room  I If  the  volume  of  the  lower  [upper]  layer  is  non-zero  (if  the  lower  [upper]  layer  volume  Is 
zero  then  the  c^^,  [Cy  ,^,]  value  is  not  used  in  the  calculation),  [(unit  of  product  K)/(kg  of  layer)] 

if  space  I is  an  'outside  room':  is  the  uniform  concentration  throughout  the  space  of  product 

of  combustion  K;  Cy  is  specified  as  being  identical  to  0^^,.  [(unit  of  product  K)/(kg  of  local 
atmosphere)] 
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^U02.l  [^U,02.ll'  I — 1 Of  2 


If  space  I is  an  inside  room:  Concentration  of  oxygen  in  lower  [upper]  layer  of  room  I if  the 
volume  of  the  lower  [upper]  layer  is  non-zero  (if  the  lower  [upper]  layer  volume  is  zero  then  the 
0^021  [0uo2il  calculation),  [(kg  of  oxygen)/(kg  of  layer)] 

If  space  I Is  an  'outside  room':  c^Qg.i  ‘S  the  uniform  concentration  throughout  the  space  of 
oxygen;  02 1 's  specified  as  being  identical  to  c^q2.i-  [(^9  of  oxygen)/(kg  of  local  atmosphere)] 


Cp 

specific  heat  at  constant  pressure  of  the  vent  flow.  [W  • s/(kg  • K)]  (suggest  1 0^  W • s/(kg  • K)  for  air 
as  default) 

NpMAX 

Maximum  allowed  number  of  products  of  combustion. 

^PROD 

Number  of  products  of  combustion,  including  oxygen,  being  tracked  in  the  simulation. 

Pdat 


Datum  absolute  pressure.  [Pa  = kg/(m  • s^] 

*^L,i  n"u,il»  I = 1 or  2 

If  space  I is  an  inside  room:  Absolute  temperature  of  the  lower  [upper]  layer  in  room  I if  the 
volume  of  the  lower  [upper]  layer  is  non-zero  (if  the  lower  [upper]  layer  volume  is  zero  then  the 
Fu.il  's  not  used  in  the  calculation).  [K] 

If  space  I Is  an  'outside  room':  T^,  Is  the  uniform  absolute  temperature  there,  taken  to  be  the 
temperature  at  the  reference  elevation,  yp^F.i;  Ty , is  specified  as  being  identical  to  T^.  [K] 

VcEiui  lyREF.J 

If  space  I is  an  inside  room:  Elevation  of  the  ceiling  [floor]  of  room  I above  the  datum  elevation, 
[m] 

If  space  I is  an  'outside  room':  yc^im  and  y^^p,  are  both  identical  and  equal  to  the  reference 
elevation  of  space  I above  the  datum  elevation,  i.e.,  the  specification  must  satisfy  ycEiui  = Vref.i’ 
The  latter  identity,  which  will  never  be  satisfied  for  an  inside  room,  is  a characteristic  of  the  input 
data  used  to  distinguish  an  inside  room  from  an  'outside  room.'  [m] 

VLAYEai’  I = or  2 

If  space  I is  an  Inside  room:  Elevation  of  the  upper/lower  layer  interface  in  room  I above  the 
datum  elevation,  [m] 

If  space  I is  an  'outside  room':  yLAVEai  specified  as  being  identical  to  yppp,.  [m] 
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VVENT 


Elevation  of  the  vent  above  the  datum  elevation.  Note  that  must  be  identical  to  either  Yceili 
or  Vrep,  for  each  of  the  one  or  two  Inside  rooms  Involved  in  the  calculation,  [m] 

^Pref.i*  J = "I  or  2 

Pressure  at  the  reference  elevation,  yR^p,,  in  space  I above  the  datum  absolute  pressure,  pp^p 
If  space  I is  an  inside  room,  then  6Prepi  and  yR^p,  must  correspond  to  the  pressure  and  elevation, 
respectively,  within  the  room  and  at  the  floor.  [Pa  = kg/(m  • s^] 


Error  tolerance  for  SpR^p,.  If  Perrorj  is  defined  as  the  uncertainty  in  6pR£p,,  I = 1 or  2,  then 
Perror.1  satisfies 

IPerror.iI  SPo  l^pREF.ll^p 

where  Po  = 1 -OPa.  The  first  term  is  based  on  an  absolute  error  tolerance  and  dominates  the 
above  error  bound  when  |6pR^p,|  is  less  than  1.0Pa.  The  second  term  is  a relative  error 
tolerance  and  dominates  when  |6pR^p,|  is  greater  than  1.0Pa.  should  be  chosen  to  be 
consistent  with  the  tolerance  specified  for  the  computation  of  ^pR^p,  terms  in  the  overall 
compartment  fire  model  computer  code  which  uses  this  algorithm. 

Pl,\  I = "I  or  2 

If  space  I is  an  inside  room:  Density  of  the  lower  [upper]  layer  in  room  I if  the  volume  of  the  lower 
[upper]  layer  is  non-zero  (if  the  lower  [upper]  layer  volume  is  zero  then  the  [p^ ,]  value  is  not 
used  in  the  calculation),  [kg/m^] 

If  space  I is  an  "outside  room":  is  the  uniform  density  there;  py , is  specified  as  being  identical 

to  Pm-  [kg/tn^] 


CALCULATION 

1.  Calculate  p,  for  I = 1 and  2 , and  then  Ap  (the  p,  calculation  follows  the  DELP  algo- 
rithm/subroutine of  reference  [8]): 


Pi  ” ^pREF.i  + ^Pi  + Pdat 

^P  ” P2  ■ Pl  “ (^PrEF,2  ■ ^pREF.l)  (^P2  " ^Pl) 


(1) 


(2) 


where 


^ yREF.I  - VvENT  - yiAYEai  ^P| PuiS^yVENT  " yREF.l) 

else:  6P|  = - PL.|9(yLAYEai  ' yREF.l)  " Pu,|9(yvENT  ' yLAYEai) 


(3) 


and  where  g,  the  acceleration  of  gravity,  is  9.8  m/s^. 
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2.  Calculate  /?,,  T,,  Cog,,.  and  the  c^,,  K = 2 to  Mppop,  for  I = 1 and  2,  and  then  Ap: 

{[(VvENT  “ yREF.l)  (VlAYER.I  “ VrEF.i)]  [(VvENT  “ yCEIUI^  (y LAYER.  I ycEiui)l}^hen: 

p\  ~ Pu.l'  '*'l  “ "^U.l'  ^02.l  “ ^U,02,l'  ^K,l  “ ^U,K,I  W 

else:  p,  = p^,,  T,  = Tm,  Cq2  , = c^Qg  ,,  and  c,^|  = C|^,^| 


Ap  = pi  - P2  (5) 

3.  Define  Vsjhigh  volume  rate  of  flow  through  the  vent,  from  the  high-  to  the  low-pressure 
space,  that  Is  predicted  with  a 'standard,'  unl-directlonal-flow-type  calculation  (i.e.,  without  regard 
to  the  effect  of  buoyancy,  in  general,  or  the  stability  of  the  cross-vent  density  configuration,  in 
particular),  where  arbitrarily  high  cross-vent  pressures  are  allowed.  Here,  the  calculation  follows 
the  model  of  [5]  as  implemented  in  the  VENTHP  algorithm/subroutine  of  [6].  Calculate  Vsjhigh- 

If  Ap  = 0 then:  Vsjhigh  = ^ CALCULATION  (6) 

If  Ap  > 0 or  Ap  < 0 then: 


a.  Define  and  compute  phigh» 

If  Ap  > 0,  i.e.,  'standard'  flow  from  (lower)  space  2 to  (upper)  space  1 , then: 

^HIGH  ~ P2»  ^ ~ ^P/P2 

If  Ap  < 0,  i.e.,  'standard'  flow  from  (upper)  space  1 to  (lower)  space  2,  then: 


^HIGH  “Pi'  ^ ^P/Pi 

(8) 

X = 1 - £ 

(9) 

Compute  C(x),  the  vent  flow  coefficient,  and  w(x): 

C(x)  = 0.85  - 0.25X  = 0.60  + 0.25e 

(10) 

w(x)  = ^ 

' 1 - [3/(47)]e  if  0 < e < 10  ® 

.f(x)/[2e]’^  if  1 >e  > 10  = 

(11) 

where 

f(x)=  ^ 

' {(27/(7  - 1)]X^[1  - if  e < 1 - [2/(7  + 

. {7[2/(7  + 1)]^  if  e > 1 - [2/(7  + 1)f^ 

(12) 

and  where  7,  the  ratio  of  specific  heats  of  the  vent  flow  gas,  is  taken  to  be  that  of  air, 
1.40.  Note  that  for  the  present  horizontal  vent  application,  the  C(x)  of  Eq.  (10)  is  taken 
to  be  consistent  with  the  standard  incompressible  limit  for  flow  through  circular  sharp- 
edged  orifices  In  the  sense  that  C -*  0.60  as  e 0 [9]. 
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(13) 


c.  Define  and  compute  Ap^^: 

APctrr  ^ [£pMAX(1  .OPa,  I^Prep^Ii  I^Pref,2I)1  ^ 

d.  Define  and  compute  F^qis^,  a numerical  damping  factor,  and  then  Vg-pHiGH* 

Fnoise  = 1.0-exp(-|Ap|’'"/Api^  (14) 

'^ST.HIGH  “ ^NOISE^W'^W(^^HIGh)  ^^|Ap|^^  (15) 

The  term  F^q,se  of  Eq.  (14)  is  designed  to  damp  out  the  numerical  noise  (error)  in  the 
calculated  value  for  Ap  that  would  othenvise  be  dominant  in  Eq.  (1 5)  when  Ap  is  small 
relative  to  the  maximum  of  1.0  pascal  and  the  calculated  reference  pressures,  6pp^p.,, 
6pREP2-  The  term  Ap^uy  of  Eq.  (13)  is  an  estimate  of  how  small  the  maximum  of  |Ap| 
must  be  to  retain  a few  digits  of  accuracy  in  the  calculation  of  Ap.  When  the  calculated 
value  of  I Ap|  Is  smaller  than  Pguj,  this  value  and,  therefore,  the  value  of  Vg-p^iQR  in  Eq. 
(15)  will  likely  contain  noise  which  should  be  damped.  F^^Qig^  is  constructed  to  tend 
towards  1 when  | Ap|  is  large  relative  to  Ap^uy  and  tends  towards  0 when  | Ap|  is  small 
relative  to  Ap^^y. 

When  the  cross-vent  density  configuration  is  unstable  (Ap  > 0),  mixed  pressure-  and  buoyancy- 
driven  aspects  of  the  flow  have  to  be  considered  and  the  vent  flow  rates  are  calculated  according 
to  [2].  Define  Vg  riqr,  Vg  low  ^ reference-[2],  buoyancy-affected,  volume  flow  rates  across 
the  vent  from  the  hlgh-to^low  and  low-to-high  pressure  spaces,  respectively. 


If  Ap  < 0:  Vg  RiQR  = Vg  LOW  = skip  to  (the  next)  item  5 of  the  CALCULATION 

a.  Calculate  T,  p,  fi(f)  in  m^/s  from  [10],  and  e > 0 : 

T = (T,  + T2)/2;  p = (p^  -h  p^/2\ 

/i(T)  = p [0.041 28(1 0'^)T®^/(T  + 110.4)];  £ = Ap/^  > 0 

b.  Calculate  D and  Gr;  if  P2  > p.,,  i.e.,  Ap  > 0,  then  replace  e by  - e < 0: 

D = Gr  = 2gD^|e|/[/i(T)/p]^;  if  Ap  > 0 then  e = - e 


(17) 


(18) 


(19) 


c.  Calculate  conditions  at  the  limit  of  uni-directional  flow  (i.e.,  the  flooding  condition)  and  the 
relative  pressure,  6p*: 

^HiGH.FLooD  ~ 0.1754Ay(2gD je|)  ^exp(0.5536e); 

^Pflood  = 0.2427(4gApD)(1  + £/2)exp(1.1072£);  (20) 

6p*  = |Ap|/AppLooD 

d.  Calculate  02  and  then  Vg  low*  Vg  h,gh: 

= Fno,seC(x)w(x)/0.1780;  02  = 1.045  (21) 
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If  5p*  > 1,  expect  uni-directional  flow: 


^B.LOW  “ ® 


~ ^HIGH.FLOODI 


,{1  - o/  + [o/  + c,=(5p*  - 1)]’'^} 
If  6p*  < 1 , expect  mixed  flow: 


(22) 


'^EX.MAX  = 0.055(4/7r)Av(gD|e|)'® 
m3  = - 0.7070;  M = (CTi/a^)^  - 1 

'^B.LOw  = Vmax[(1  + ni3/2)(1  - <5P*)^  - (2  + mJ2W  - 5p*)]= 

— '^B.LOW  + {M  ' [1  + - 1)(1  - ^P*)l'^}'i*HIGH.FLOOc/(^  " 


(23) 

(24) 


(Note:  of  Eq.  (21)  Is  somewhat  different  than  in  [2].  The  modification  here  is 

provide  for  an  analytic  representation  of  vent  flow  which  is  continuous  and  uniformly  valid 
even  as  6p’'  is  increased  to  a level  where  compressibility  effects  become  important.  The 
result  of  [2]  does  not  include  the  effect  of  compressibility.) 


Following  [2]^  when  Lp  > 0 the  above  will  be  taken  as  the  flow  solution  provided  Gr  > 2(10^). 
For  smaller  Gr,  the  reference-[2]  sojution  begins  to  loose  its  validity,  and  there  is  no  existing 
model  for  the  flows  the  range  0 < Gr  <_2(10^).  However,  it  is  clear  that  the  above-calculated 
■standard'  flow  must  be  approached  as  Gr  0. 


Consistent  with  the  above  comments,  when  Gr  is  in  the  range  0 < Gr  < 2(10^),  the  vent  flow  will 
be  estimated  by: 

flow  = ('standard'  flow)  + [(reference-[2]  flow)  - ('standard'  flow)]Gr/[2(10^)] 

a Define  and  Vlqw  ^ volume  flow  rates  across  the  vent  from  the  high-to-low  and 
low-to-high  pressure  spaces,  respectively. 

If  Ap  < 0:  = VgjHiGH*  ^Low  = ® (25) 

If  Ap  > 0 and  0 < Gr  < 2(10^): 

^HIGH  ” ^ST.HIGH  (^B.HIGH  " ^ST.HIGh)®‘^^I2(‘^  0^)1 

(26) 

^Low  “ ^b.low®*^^[2(1  0 )] 

If  Ap  > 0 and  Gr  > 2(10^):  Vhiqh  = ^b.highJ  ^low  “ ^b,low  (27) 

b.  Define  and  calculate  I = 1 and  2,  as  the  vent  volume  flow  rate  entering  space  I. 

If  Ap  > 0:  = V^ighI  ^vent,2  ” ^low  (28) 

(i.e.,  flow  to  the  upper  space  is  the  flow  from  the  high-  to  the  low-pressure  space,  and 
flow  to  the  lower  space  Is  the  flow  from  the  low-  to  the  high-pressure  space) 
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If  Ap  < 0:  VvENTT.I  “ ^LOW>  ^VENT2  ~ ^HIGH  (29) 

(i.e.,  flow  to  the  upper  space  is  the  flow  from  the  low-  to  the  high-pressure  space,  and 
flow  to  the  lower  space  is  the  flow  from  the  high-  to  the  low-pressure  space) 

Calculate  the  vent  flow  properties  "^vent.i’  ^vent.02,i’  ^vent,k,i’  K ~ 2 to  NppQQ,  I = 1 and  2! 

^VENT.I  “ P2(Pl/P2)>  ^VENT.2  ~ "^VENT.I  ~ "^2’  '^VENT,2  ~ (^O) 

^ENT,02.1  ~ ^02,2'  ^Errr.02.2  “ ^02,1  > ^ENT.K.1  ” ^K,2'  ^ENT,K.2  ~ ^K.1>  K = 2 tO  NppQp 

Calculate  the  vent  flow  rates  I = 1 and  2: 

^VENTJ  “ PVENT.rVENT.I  (21) 

Calculate  the  vent  flow  rates  ^02venti’  ^ kventi*  ^ ” 2 to  NppQQ,  I ~ 1 and  2. 

^VENT.I  “ ^VENT.I^p'^VENT.P  ^02, VENT.  I = ^VENTI^ENT.02.P 

(32) 

^K.VENTI  “ ^VENTI^ENT.K.I’  K = 2 tO  NppQp 

Calculate  rates  at  which  flows  are  added  to  layers  of  each  space  as  a result  of  the  vent  flow 
extracted  from  It.  First  consider  space  1 (I  = 1,  J = 2)  and  then  space  2 (I  = 2,  J = 1).  For  either 
case: 


{[(VvENT  ” VrEF.i)  (VlAYEai  “ VrEF.i)]  O'"  [(VvENT  ” VcElUl)  (VlAYER,!  yCEIL,|)]}  (‘•®-» 

the  vent  flow  to  room  J is  extracted  from  the  upper  layer  of  room  I and,  if  space  I is  an 
inside  room,  the  lower  layer  of  room  I Is  unchanged)  then: 

My  I = - ^ui  ” ^u.i  “ ’ ^vent.J’  ^l.i  “ 

(33) 

^02.U.I  ~ ■ ^02.VENTJ'  ^02.LI  “ ^K.U.1  “ " ^KVENT.J’  ^K.UI  “ ® 

If  y^gp,  = (i.e.,  space  I is  an  outside  space)  modify  the  results  of  Eqs.  (33) 

as  follows: 


M 


ui 


M 


u.i» 


^ui  ” ^u.i' 


p sz  p ‘ p =1  p 

^02.LI  ^02.U.I'  ^K,L,I  ^K.U.1 


(34) 


If  the  condition  above  Eqs.  (33)  is  not  satisfied  (i.e.,  the  vent  flow  to  room  J is  extracted 
from  the  lower  layer  of  room  I and,  if  space  I is  an  inside  room,  the  upper  layer  of  room 
I is  unchanged)  then: 


^L,l ^VENT.J’  ^u.l  “ ^ui  “ ■ ^VENT.J’  ^U.l  “ 

(35) 

^02.1,1  “ ■ ^02.VENT,J'  ^02,U.I  “ 9;  P,^L.|  = - Pk^veNT.J'  ^K.U.1  ” 9 

If  ypppi  = ycEiui  (i*®M  space  I is  an  outside  space)  modify  the  results  of  Eqs.  (35) 
as  follows: 


Mu.,  = M 


L,l* 


‘U.l 


*L.I'  '^02.U.I 


= P 


02.UI*  '^K.U.I 


= P 


K.UI 


(36) 
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SUBROUTINE  VARIABLES 

All  nomenclature  in  the  subroutine  is  identical  to  the  nomenclature  used  above  except  for: 
Ay  - AVENT  [m^j 

C - COEF 
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Cp 

CP  [W.s/(kg.K)] 

^LK.I'  ^U,K.I 

C0NL(K,l),  C0NU(K,l),  1 = 1 or  2 [(unit  of  product  K)/(kg 
of  layer)] 

^L.02,1'  ^U,02.l 

C0NL(1,I),  C0NU(1,I),  I = 1 or  2 [(kg  of  oxygen)/(kg  of 
layer)] 

^ENT,K.I 

CVENT(K,I),  l=  1 or  2 [(unit  of  product  K)/(kg  of  vent 
flow) 

^ENT.02.1 

CVENT(1,I),  I = 1 or  2 [(kg  of  oxygen)/(kg  of  vent  flow)] 

Ck.1 

C(K,I),  l=  1 or  2 [(unit  of  product  K)/(kg  of  vent  flow) 

^02.1 

C(1,l),  I = 1 or  2 [(kg  of  oxygen)/(kg  of  vent  flow)] 

^NOISE 

FNOISE  [dimensionless] 

f 

FF  [dimensionless] 

Gr 

GR  [dimensionless] 

g 

9.8  m/s^ 

M 

XM  [dimensionless] 

Mj , 

XML(I),  XMU(I),  I = 1 or  2 [kg/s] 

^VENT.I 

XMVENT(I).  1 = 1 or  2 [kg/s] 

rtig 

XM3  [dimensionless] 

*^piyiAx 

NPMAX  [dimensionless] 

^PROD 

NPROD  [dimensionless] 

^KUI’  ^KU.I 

PL(K,I),  PU(K,I),  1 = 1 or  2 [(unit  of  product  K)/s] 

b 

^K.VENT.1 

PVENT(K,I),  1 = 1 or  2 [(unit  of  product  K)/s] 

b b 

^02.UI»  ^02.U.I 

PL(1,I),  PU(1,I),  I = 1 or  2 [(kg  of  oxygen)/s] 

b 

^02.VEMT.I 

PVENT(1,I),  1 = 1 or  2 [(kg  of  oxygen)/s] 

Pdat 

PDATUM  [Pa  = kg/(m  s^] 

Pi 

P(l),  1 = 1 or  2 [Pa  = kg/(m-s^] 

O 

r 

O 

c 

QL(I),  QU(I).  1 = 1 or  2 [W] 

^VENT.I 

QVENT(I).  1 = 1 or  2 [W] 
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T 

TBAR  [K] 

T, 

T(l),  1 = 1 or  2 [K] 

\l>  "^U.l 

TL(I),  TU(I),  I = 1 or  2 [K] 

*^VENT,l 

TVENT(I),  I = 1 or  2 [K] 

^B.HIGH 

VBHIGH  = VB(1)  if  Ap  > 0,  = VB(2)  if  Ap  < 0 [m=^/s] 

^B.LOW 

VBLOW  = VB(2)  if  Ap  > 0,  = VB(1)  if  Ap  < 0 [m^/s] 

^EX,MAX 

VEXMAX  [m^/s] 

^LOW 

VHIGH,  VLOW  [m^/s] 

^HIGH.FLOOD 

VHIGHFL  [m^/s] 

^ST.HIGH 

V = VST(1)  if  Ap  > 0,  = VST(2)  if  Ap  < 0 [m^/s] 

^VENT.I 

WENT(I),  1 = 1 or  2 [m^/s] 

W 

W [dimensionless] 

X 

X [dimensionless] 

VCEIUI 

YCEIL(I),  1 = 1 or  2 [m] 

yLAYER.1 

YLAY(I),  1 = 1 or  2 [m] 

Vref.i 

YREF(I),  1 = 1 or  2 [m] 

VvENTT 

YVENT  [m] 

7 

1.40 

Ap 

DELP  [Pa  = kg/(m.s^] 

DPC1 D2  [Pa  = kg/(m  • s^] 

^Pflood 

DELPFD  [Pa  = kg/(m-s^] 

Ap 

DELDEN  [kg/m^] 

5p* 

DPDDPFL  [Pa  = kg/(m  • s^] 

6P| 

DP(I),  1 = 1 or  2 [Pa  = kg/(m.s^] 

^Pref.i 

DPREF(I),  1 = 1 or  2 [Pa  = kg/(m-s^] 

e 

EPS  [dimensionless] 
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^HIGH 

P 

P\ 

Pl,\'  Pu,\ 
Advent,  I 

PREPARED  BY 

Leonard  Y Cooper 
December  1993 


EPSP  [dimensionless] 

XMEW  [m^/s] 

DENHIGH  [kg/m^] 

DENBAR  (kg/m®l 
DEN(I),  I = 1 or  2 [kg/m^] 

DENL(I),  DENU(I),  I = 1 or  2 [kg/m^] 
DENVNT(I),  I = 1 or  2 [kg/m'*] 
SIGMA1  [dimensionless] 

SIGMA2  [dimensionless] 
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(C) 


Figure  1 . The  possible  configurations  of  the  two  spaces  Joined  by  a horizontal  ceiling/floor  vent 
with  space  1 above  space  2:  a)  two  inside  rooms;  b)  an  outside  space  above  an 
inside  room;  c)  an  inside  room  over  an  outside  space. 
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Figure  2.  The  geometry  and  conditions  local  to  a horizontal  celling/floor  vent  which  determine 
the  characteristics  of  the  vent  flow. 


VENTCF2 - 16 


oooooooooooooooooooooooooooooooooooooooooooooooooooo 


SUBROUTINE  VENTCF2 


SUBROUTINE  VENTCF2( 

I AVENT,CONL,CONU,CP,NPMAX,NPROD,PDATUM,TL,TU, 

I YCEIL,YREF,YLAY,YVENT,DPREF,EPSP,DENL,DENU, 

O C,CVENT^ML,XMU,XMVENT,PL,PU,PVENT,P,QL,QU,QVENT,T,TVENT, 

O DELP.DEN.DENVNT) 

INCLUDE  PRECIS.INC 
C*BEQ 

C***  VENTCF  - CALCULATION  OF  THE  FLOW  OF  MASS,  ENTHALPY,  OXYGEN  AND 
C OTHER  PRODUCTS  OF  COMBUSTION  THROUGH  A HORIZONTAL 

C VENT  JOINING  AN  UPPER  SPACE  1 TO  A LOWER  SPACE  2.  THE 

C SUBROUTINE  USES  INPUT  DATA  DESCRIBING  THE  TWO-LAYER 

C ENVIRONMENT  OF  INSIDE  ROOMS  AND  THE  UNIFORM 

C ENVIRONMENT  IN  OUTSIDE  SPACES. 

C 

c*«*  SUBROUTINE  ARGUMENTS 


INPUT 


AVENT 

CONL(K,l) 


CONL(K,0 


CP 

NPMAX 

NPROD 

PDATUM 

TL(0 

TU(0 

YCEILO) 

YREF(0 

YLAY(I) 

WENT 

DPREFO) 

EPSP 

DENL(0 

DENU(I) 


OUTPUT 


AREA  OF  THE  VENT  [M**21 

CONCENTRATION  OF  PRODUCT  K IN  LOWER  LAYER  OF  AN  INSIDE 
ROOM  I OR  UNIFORM  CONCENTRATION  OF  PRODUCT  K IN  AN 
OUTSIDE  SPACE  I [(UNIT  OF  PRODUCT)/(KO  LAYER)] 
CONCENTRATION  OF  PRODUCT  K IN  UPPER  LAYER  OF  AN  INSIDE 
ROOM  I OR  UNIFORM  CONCENTRATION  OF  PRODUCT  K IN  AN 
OUTSIDE  SPACE  I [(UNIT  OF  PRODUCT)/(KO  LAYER)] 

SPEaFIC  HEAT  [W*S/(KO*K)] 

MAXIMUM  ALLOWED  NUMBER  OF  PRODUCTS 
NUMBER  OF  PRODUCTS  IN  CURRENT  SCENARIO 
DATUM  ABSOLUTE  PRESSURE  [PA  > KO/(M*S**2)] 

TEMPERATURE  OF  LOWER  LAYER  OF  AN  INSIDE  ROOM  I OR 
TEMPERATURE  OF  AN  OUTSIDE  SPACE  I [K] 

TEMPERATURE  OF  UPPER  LAYER  OF  AN  INSIDE  ROOM  I OR 
TEMPERATURE  OF  AN  OUTSIDE  SPACE  I [K] 

HEIGHT  OF  CEIUNG  ABOVE  DATUM  ELEVATION  FOR  AN  INSIDE 
ROOM  I OR  YREF(0  FOR  AN  OUTSIDE  SPACE  I [M] 

HEIGHT  OF  REFERENCE  ELEVATION  FOR  SPACE  I ABOVE  DATUM 
ELEVATION  [M] 

HEIGHT  OF  LAYER  ABOVE  DATUM  ELEVATION  FOR  AN  INSIDE 
ROOM  I OR  YREF(0  FOR  AN  OUTSIDE  SPACE  I [M] 

HEIGHT  OF  VENT  ABOVE  DATUM  ELEVATION  [M] 

PRESSURE  IN  SPACE  I AT  ITS  REFERENCE  ELEVATION  ABOVE 
DATUM  ABSOLUTE  PRESSURE  [PA  « KO/(M*S**2)] 

ERROR  TOLERANCE  FOR  DPREF  [DIMENSIONLESS] 

DENSITY  OF  LOWER  LAYER  OF  AN  INSIDE  ROOM  I OR  DENSITY 
OF  AN  OUTSIDE  SPACE  I [KO/M**3] 

DENSITY  OF  UPPER  LAYER  OF  AN  INSIDE  ROOM  I OR  DENSITY 
OF  AN  OUTSIDE  SPACE  I [KO/M**3] 


C(K,0 

CVENTO) 

XML(I) 

XMU(I) 

XMVENT(I) 

PL(K,I) 


CONCENTRATION  OF  EACH  PRODUCT  IMMEDIATELY  ABOVE  (IN 
SPACE  I B 1)  AND  BELOW  (IN  SPACE  I = 2)  THE  VENT 
[(UNIT  OF  PRODUCT)/(KO  LAYER)] 

CONCENTRATION  OF  EACH  PRODUCT  IN  THE  VENT  FLOW 
COMPONENT  ENTERING  SPACE  I [(UNIT  OF  PRODUCT)/(KO  OF 
VENT  FLOW)] 

RATE  AT  WHICH  MASS  IS  ADDED  TO  THE  LOWER  LAYER  OF  AN 
INSIDE  ROOM  I - 1 (I  > 2)  OR  TO  AN  OUTSIDE  SPACE  I * 

1 (I  > 2)  DUE  TO  THE  VENT  FLOW  COMPONENT  ENTERING 
SPACE  « 2 (I  > 1)  [KG/S] 

RATE  AT  WHICH  MASS  IS  ADDED  TO  THE  UPPER  LAYER  OF  AN 
INSIDE  ROOM  I - 1 0 - 2)  OR  TO  AN  OUTSIDE  SPACE  I - 1 
(1-2)  DUE  TO  THE  VENT  FLOW  COMPONENT  ENTERING  SPACE 
1-20-1)  [KG/S] 

MASS  FLOW  RATE  IN  THE  VENT  FLOW  COMPONENT  ENTERING 
SPACE  I [KG/S] 

RATE  AT  WHICH  PRODUCT  K 18  ADDED  TO  THE  LOWER  LAYER  OF 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C*END 


PU(K,I) 

PVENT(K,i) 

PO) 

QL(1) 

QUO) 

OVENTO) 

TO) 

TVENTO) 

DELP 

DENO) 

DENVNTO) 


AN  INSIDE  RCX)M  I = 1 (I  2)  OR  TO  AN  OUTSIDE  SPACE  I 
= 1 (I  = 2)  DUE  TO  THE  VENT  FLjOW  CONPONENT  ENTERING 
SPACE  I s 2 (I  = 1)  [(UNIT  OF  PRODUOO/S] 

RATE  AT  WHICH  PRODUCT  K IS  ADDED  TO  THE  UPPER  LAYER  OF 
AN  INSIDE  ROOM  I = 1 0 » 2)  OR  TO  AN  OUTSIDE  SPACE  I 
« 1 0 = 2)  DUE  TO  THE  VENT  FLOW  COMPONENT  ENTERING 
SPACE  I = 2 (I  = 1)  [(UNIT  OF  PRODUCT)/S] 

FLOW  RATE  OF  PRODUCT  K IN  THE  VENT  FLOW  COMPONENT 
ENTERING  SPACE  I [(UNIT  OF  PROOUCT)/S] 

ABSOLUTE  PRESSURE  IMMEDIATELY  ABOVE  (IN  SPACE  I > 1) 

AND  BELOW  ON  SPACE  I ^ 2)  THE  VENT  [PA  ^ KG/(M*S**2)] 

RATE  AT  WHICH  ENTHALPY  IS  ADDED  TO  THE  LOWER 
LAYER  OF  AN  INSIDE  ROOM  I = 1 0 « 2)  OR  TO  AN  OUTSIDE 
SPACE  I c 1 0 *=  2)  DUE  TO  THE  VENT  FLOW  COMPONENT 
ENTERING  SPACE  I 2 0 = 1)  [W] 

RATE  AT  WHICH  ENTHALPY  IS  ADDED  TO  THE  UPPER  LAYER  OF 
AN  INSIDE  ROOM  I = 1 (I  « 2)  OR  TO  AN  OUTSIDE  SPACE  I 
« 1 0 » 2)  DUE  TO  THE  VENT  FLOW  COMPONENT  ENTERING 
SPACE  I = 2 (I  B 1)  [W] 

FLOW  RATE  OF  ENTHALPY  IN  THE  VENT  FLOW  COMPONENT 
ENTERING  SPACE  I [W] 

ABSOLUTE  TEMPERATURE  IMMEDIATELY  ABOVE  ON  SPACE  I « 

1)  AND  BELOW  (IN  SPACE  I « 2)  THE  VENT  [K] 

ABSOLUTE  TEMPERATURE  OF  THE  VENT  FLOW  COMPONENT 
ENTERING  SPACE  I [K] 

CROSS-VENT  PRESSURE  DIFFERENCE,  P(2)  - P(1),  [PA  « 
KG/(M*S**2)1 

DENSITY  IMMEDIATELY  ABOVE  ON  SPACE  I = 1)  AND  BELOW 
ON  SPACE  I « 2)  THE  VENT  [KO/M**3] 

DENSITY  OF  THE  VENT  FLOW  COMPONENT  ENTERING  SPACE  I 
[KG/M**31 


C****  NOTE  THAT  NPMAX2  SHOULD  BE  BIGGER  THAT  NPMAX 
PARAMETER  (NPMAX2-10) 

DIMENSION  CONL(NPMAX2,2),  CONU(NPMAX2,2),  TL(2),  TU(2) 

DIMENSION  YCEIL(2),  YREF(2),  YLAY(2),  DPREF(2) 

DIMENSION  DENL(2).  DENU(2) 

DIMENSION  C(NPMAX2,2),  CVENT(NPMAX2,2),  XML(2),  XMU(2) 

DIMENSION  XMVENT(2) 

DIMENSION  PL(NPMAX2,2).  PU(NPMAX2,2),  PVENT(NPnUO(2,2).  P(2) 

DIMENSION  QL(2),QU(2) 

DIMENSION  QVENT(2),  T(2).  TVENT(2).  DEN(2),  DENVNT(2) 

DIMENSION  DP(2),VB(2),VST(2),  VVENT(2) 

PARAMETER  (QAM*i1.40D0) 

DATA  IHRST/O/ 

SAVE  IFIRST.QAMCUT,OAMMAX 
PARAMETER  (O>9.80D0) 

MRAMETER  (PI«3.1415926S4D0) 

MRAMETER  (SIGMA2-1.04SD0) 

PARAMETER  (XM3s-O.7070D0) 

QR«0.00 
VHIQHFL«O.DO 
DELPFL-O.DO 
DO  SI-1,2 
VST0)*O.DO 
VB(l)-O.DO 
OVENTO)  >0.D0 
XMVENT0)«O.DO 
DO  2 NP-1.NPROD 
PVENT(NP,l)-O.DO 
2 CONTINUE 

5 CONTINUE 

C 

C***  THE  FOLLOWING  CODE  SEGMENT  COMPUTES  CONSTANTS  REQUIRED  BY  VENTCF2. 
C***  IT  IS  EXECUTED  THE  FIRST  TIME  VENTCF  IS  CALLED. 

C 
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c 

C***  1. 

c*** 

c 


c 

C*** 

c*** 

c*** 

c 


c 

c*** 

c*** 

c*** 

c 

$ 


10 

c 

c*** 

c*** 

c 


22 

$ 

$ 


24 


26 

30 

C 

c*** 

c*** 

c 

c 

c*** 

c 

c*** 

c 


IF(IFIRST.EQ.0)THEN 
FIRST  * 1 

OAMCUT  s (2.0D0/(QAM+1.0D0))**(QAM/(QAM-1.0D0)) 
ZZZ=QAM*((2.0D0/(QAM +1 .0D0))**((QAM  -«-1 .0D0)/(QAM-1  .ODO))) 
OAMMAX  « DSQRT(ZZZ) 

ENDIF 

AND  2.  CALCULATE  THE  P(l),  DELP,  THE  OTHER  PROPERTIES 

ADJACENT  TO  THE  TWO  SIDES  OF  THE  VENT,  AND  DELDEN. 

DO  10  I - 1,  2 

IF(YREFO).LE.YVENTAND.YVENT.LE.YLAY(l))THEN 

THE  VENT  IS  AT  OR  BELOW  THE  REFERENCE  ELEVATION  IN  SPACE  I.  IF 
SPACE  I IS  AN  INSIDE  ROOM  THEN  BOTH  THE  VENT  AND  THE  LAYER 
INTERFACE  ARE  AT  THE  FLOOR  ELEVATION. 

DP(I)  * - 0*DENL(l)*(YVENT  - YREF(I)) 

ELSE 

THE  VENT  IS  ABOVE  THE  REFERENCE  ELEVATION  IN  SPACE  I. 

IF  SPACE  I IS  AN  INSIDE  ROOM  THEN  THE  VENT  IS  AT  THE 
CEIUNQ. 


DP(I)  = - G*DENL(0*(YLAY(l)-YREF(l)) 

- 0*DENU(l)*(YVENT-YLAY(l)) 

ENDIF 

P(0  « DPREF(I)  + DP(I)  + PDATUM 
CONTINUE 

DELP  IS  PRESSURE  IMMEDIATELY  BELOW  THE  VENT  LESS  PRESSURE 
IMMEDIATELY  ABOVE  THE  VENT. 

DELP  « (DPREF(2)-DPREF(1))+PP(2)-DP(1)) 

DO  30  I s 1.2 

DEN(0  E O.ODO 

T(l)  > 0.0D0 

DO  22  K « 1.  NPROD 

C(K,0  > O.ODO 
CONTINUE 

IF(((YVENT.EQ.YREF(0)AND.((YLAY(0-YREF(I)).LT.O.DO)) 

.OR.((YVENT.EQ.YCEIL(I))AND.((YCEIL(0-YLAY(I)) 

.QT.0.D0)))THEN 

DEN(0  « DENU(0 

T(I)  -TU{0 

DO  24  K s 1.  NPROD 

C(K,I)  > CONU(K.O 
CONTINUE 

ELSE 

DEN(Q  « DENLO) 

T(I)  « TL(0 

DO  26  K « 1,  NPROD 

C(K.O  - CONL(K.O 
CONTINUE 

ENDIF 

CONTINUE 

DELDEN  IS  DENSITY  IMMEDIATELY  ABOVE  THE  VENT  LESS  DENSITY 
DENSITY  IMMEDIATELY  BELOW  THE  VENT 

DELDEN  - DEN(1)  • DEN(2) 

3.  CALCULATE  VST(0.  THE  "STANDARD"  VOLUME  RATE  OF  FLOW 
THROUGH  THE  VENT  INTO  SPACE  I 

CALCULATE  VST(0  IF  DELP  « 0 

IF(DELP.Ea0.0D0)THEN 
VST(1)  « 0.0D0 
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UULUo  oUULi 


VST(2)  »=  O.ODO 

ENOIF 

IFPELP.EQ.O.DO)  GOTO  32 
C 

C***  CALCULATE  VST(I)  FOR  NONZERO  DELP 
C 

IFPELP.QT.0.0D0)THEN 

VST(2)  = 0.0D0 
RHO  = DEN(2) 

EPS  = DELP/P(2) 

ENOIF 

IFPELP.LT.0.0D0)THEN 

VST(1)  s O.ODO 
RHO  » DEN(1) 

EPS  « -DELP/PO) 

ENOIF 

X s 1.000  - EPS 

COEF  = 0.6000  + 0.2SD0*EPS 

XO  « 1.000 

EPSCUT  = EPSP*MAX(XO,DPREF(1),DPREF(2)) 

EPSCUT  « DSQRT(EPSCUT) 

SROELP  s DSQRT(DABSPELP)) 

FNOISE  s 1.000 

IF((SRDELP/EPSCUT).LE.130.D0)THEN 

FNOISE  nr  1.000  • DEXP(-SRDELP/EPSCUT) 

C 

C***  NOTE:  IF  SINGLE  PRECISION  THEN  USED  65.  INSTEAD  OF  130. 

C 

ENOIF 

IF(EPS.LE.0.1  D-5)THEN 

W ar  1.000  • 0.7SD0*EPS/OAM 

ELSE 

IF(EPS.LT.OAMCUT)THEN 

GO  « X**(1.0D0/GAM) 

FF  E DSQRT((2.0D0*GAM/(GAM-1.0D0))*GG*GG* 

$ (1.0DO-X/GG)) 

ELSE 

FF  « GAMMAX 

ENOIF 

W » FF/DSQRT(EPS+EPS) 

ENOIF 

V - FNOISE*COEF*W*DSQRT(2.0O0/RHO)*AVENT*SRDELP 
IFPELP.GT.O.OOO)VST(1)  » V 
IF(DELP.LT.0.0D0)VST(2)  V 
32  CONTINUE 
C 

4.  WHEN  CROSS^NT  DENSITY  CONFIGURATION  IS  UNSTABLE,  I.E., 

DELDEN  > 0,  THEN  CALCULATE  THE  VENT  FLOW  ACCORDING  TO: 

COOPER,  LY.  COMBINED  PRESSURE-  AND  BUOYANCY-DRIVEN  FLOW 
THROUGH  A HORIZONTAL  VENT.  TO  APPEAR  AS  NISTIR,  NATIONAL 
INSTITUTE  OF  STANDARDS  AND  TECHNOLOGY  GAITHERSBURG  MO 

FOR  STABLE  CONFIGURATION.  GO  TO  5.  AND  VOLUME  FLOW  RATES  WITH 
STANDARD  MODEL 

IF(DELDEN.LE.0.0D0)  GOTO  35 

FOR  UNSTABLE  CONRQURATION,  NOW  CALCULATE  THE  COMBINED  PRESSURE- 
AND  BUOYANCY-DRIVEN  VOLUME  VENT  FLOW  RATES  FROM  THE  HIGH-TO-LOW 

AND  FROM  THE  LOW-TO-HIGH  SIDES  OF  THE  VENT,  VBHIGH  AND  VBLOW, 
RESPECTIVELY  FROM  THESE  THEN  CALCULATE  THE  COMBINED  PRESSURE- 
AND  BUOYANCVORIVEN  VOLUME  FLOW  RATES.  VB(I),  THROUGH  THE  VENT 
INTO  SPACE  I. 

C 

TBAR«  (T(1) -f  T(2))/2.D0 
DENBAR>  PEN(1)  •l-DEN(2))/2.D0 

XMEW-(0.04128D-7)*DENBAR*(TBAR**2.SD0)/(TBAR-I-110.4D0) 

EPSDEN«DELDEN/DENBAR 

D>DSQRT(4.DO*AVENT/P0 
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1 

1 


1 


1 


1 

1 


c 

Q*** 

Q*** 

c 

35 


C*** 

c 


50 

c 

c*** 

c 


60 

c 

Q*** 

c 


OR=2.DO*Q*(0**3)*EPSDEN/((XMEW/DENBAR)**2.DO) 

IFPELRQT.0.D0)  EPSOEN=-EPSDEN 

VHIQHFL=0.1754D0*AVENT*DSQRT(2.D0*O*D*DABS(EPSDEN))* 

DEXP(0.5536D0*EPSDEN) 

DELPFL=0.2427DO*(4.DO*Q*OABS(EPSDEN*DENBAR)*D)*(1.DO-i-EPSDEN/2.DO)* 

DEXP(1 .1 072D0*EPSDEN) 

DPDDPFL= DABS(DELP)/DELPFL 
SIOMA1  sFNOISE*COEF*W/0.1780D0 
iFpPODPFLQE.1  .D0)THEN 
VBLOW^O.DO 

VBHIQH  = VHIQHFL*(1  .D0-(SIQMA2*^.D0)  + 

DSQRT(SIGMA2**4.D0+ (SIOMA1**2.DO)*  pPDDPFL-1  .DO))) 

ELSE 

VEXMAX=0.0S5D0*(4.D0/PI)*AVENT*DSQRT(O*D*DABS(EPSDEN)) 

XM = (SIQMA1  /SIQMA2)**2.D0-1  .DO 

VBLDW=VEXMAX*(((1  .D0+XM3/2.D0)*((1  .D0-DPDDPFL)**2)- 

(2.D0+XM3/2.D0)*(1.D0-DPDDPFL))**2.D0) 

IF(DPDDPFLEQ.O.DO)THEN 
VBHIQH =VBLOW 

ELSE 

VBHIQHbVBLOW-I- 

(XM-DSORT(1  .D0+ (XM**2.D0-1  .D0)*(1  .DO-DPDDPFL)))* 

VHIQHFL/(XM-1.D0) 

ENDIF 

ENDIF 

IFPELP.QT.0.D0)THEN 

VB(1)sVBHIQH 

VB(2)sVBLOW 

ELSE 

VB(1)=VBLOW 
VB(2)s  VBHIQH 

ENDIF 

5.  CALCULATE  VVENT(0,  THE  VOLUME  RATE  OF  FLOW  THROUGH 
THE  VENT  INTO  SPACE  I 

CONTINUE 

IFPELOEN.LE.O.ODO)THEN 

VVENT(I)  = VST(1) 

VVENT(2)  = VST(2) 

ELSE 

IF((O.DO.LT.QR)JVND.(OR.LT.2.D7))THEN 

VVENT(1)«VST(1)+(OR/2.D7)*(VB(1)-VST(1)) 

VVENT(2)  * VST(2) + (QR/2.D7)*(VB(2)-VST(2)) 

ELSE 

VVENT(1)=VB(1) 

VVENT(2)bVB(2) 

ENDIF 

ENDIF 

6.  CALCULATE  THE  VENT  FLOW  PROPERTIES 

DENVNT(I)  - DEN(2)*P(1)/P(2) 

DENVNT(2)  « DEN(1)*P(2)/P(1) 

TVENT(I)  - T(2) 

TVENT(2)  « T(1) 

DO  50  K > 1,  NPROD 

CVENT(K,1)  * C(K^) 

CVENT{K^)  » C(K,1) 

CONTINUE 

7.  CALCULATE  THE  VENT  MASS  FLOW  RATES 
DO60 1 - i;z 

XMVENTO)  - DENVNT{0*VVENT(0 
CONTINUE 

8.  CALCULATE  THE  REST  OF  THE  VENT  FLOW  RATES 
DO  70  l>1^ 
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QVENT(I)  = XMVENT(I)*CP*TVENT(I) 

D0  65  K = 1,  NPROD 

PVENT(K.I)  « XMVENT(I)*CVENT(K.I) 

65  CONTINUE 

70  CONTINUE 

C 

C***  9.  CALCULATE  THE  RATE  AT  WHICH  THE  VENT  FLOWS  ADD  MASS, 

C***  ENTHALPY,  AND  PRODUCTS  TO  THE  LAYERS  OF  THE  SPACES  FROM 

C***  WHICH  THEY  ARE  EXTRACTED.  FIRST  TREAT  LAYERS  OF  SPACE  1 

C***  (1=1,  J=2)  AND  THEN  SPACE  2 (1=2,  J=1). 

C 

DO  95  I = 1,2 

IF  (I.EQ.1)THEN 
J = 2 

ELSE 

J = 1 

ENDIF 

IF(((YVENT.EQ.YREF0)).AND.((YLAY(l)-YREF(l)).LT.O.DO)) 

$ .OR.((YVENT.EQ.YCEIL(l))JU4D.((YCEIL(l)-YLAY(l)) 

$ .QT0.D0)))THEN 

XMU(I)  = -XMVENT(J) 

XML(I)  = 0.0D0 
QU(I)  = -aVENT(J) 

QL(I)  = 0.0D0 
DO  75  K = 1.  NPROD 

PU(K,I)  = -PVENT(K,J) 

PL(K,I)  = 0.0D0 

75  CONTINUE 

IFPABS(YREF(0-YCEIL(0).LT.0.0001  D0)THEN 
XML(I)  = XMU(I) 

QL(I)  = QU(0 
DO  80  K = 1,  NPROD 

PL(K,I)  = PU(K,0 

80  CONTINUE 

ENDIF 

ELSE 

XML(0  = •XMVENT(J) 

XMU(i)  = 0.0D0 
QLO)  = -QVENT(J) 

QU(I)  = 0.0D0 
DO  85  K = 1,  NPROD 

PL(K.I)  = •PVENT(K,J) 

PU(K,I)  = 0.0D0 

85  CONTINUE 

IF(OABS(YREF(0-YCEIL(0).LT.0.0001  D0)THEN 
XMU(0  « XML(I) 

QUO)  = QL(0 
DO  90  K = 1,  NPROD 

PU(K.I)  = PL(K.I) 

00  CONTINUE 

ENDIF 

ENDIF 

95  CONTINUE 

RETURN 
END 
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VENTCF2A  - CALCULATION  OF  THE  FLOW  THROUGH  A HORIZONTAL  CEILING/FLOOR  VENT 
CONNECTING  TWO  SPACES  WITH  "SMOOTHING"  OF  LAYER  EXTRACTION  RATES 
AT  TIMES  OF  RELATIVELY  THIN  ADJACENT-VENT  LAYERS 


DESCRIPTION 

Consider  an  instant  of  time  during  the  simulation  of  a multi-room  compartment  fire  environment.  This 
algorithm  calculates  the  flow  of  mass,  enthalpy,  oxygen,  and  other  products  of  combustion  through  a 
horizontal  vent  located  In  a ceiling/floor  partition  common  to  any  two  inside  rooms  of  the  facility  or 
between  an  inside  room  and  the  outside  environment  local  to  the  vent.  VENTCF2A  is  a modification  of 
VENTCF2  [1].  It  provides  special  considerations  for  "smoothing"  rates  of  layer  extraction  from  the  flow 
source  room  at  times  of  relatively  thin  adjacent-vent  layers.  When  used  in  a full  zone  model,  and 
depending  on  the  integration  software  for  the  particular  model,  the  considerations  in  VENTCF2A  eliminate 
singularities  that  may  cause  convergence  problems  In  fire  simulations  at  times  when  adjacent-vent  layers 
are  growing  or  shrinking  from  near-zero  depths. 

Depicted  in  Figure  1 a Is  the  vent  and  the  two  spaces  when  they  are  both  Inside  rooms  of  a multi-room 
facility.  Figures  1 b and  1 c depict  the  situation  when  the  two  spaces  Involve  one  Inside  room  of  the  facility 
and  one  outside  space,  either  above  or  below  the  room,  in  which  is  simulated  the  outside  environment 
local  to  the  vent. 

As  in  Figure  1,  designate  the  top  space  as  space  1 and  the  bottom  space  as  space  2.  It  is  assumed  that 
the  temperature,  density,  concentration  of  oxygen  and  of  other  products  of  combustion  of  Interest  In  the 
upper  and  lower  layer  of  each  inside  room  and  in  the  environment  local  to  the  vent  of  an  outside  space 
are  specified.  Also  specified  in  each  inside  room  are:  the  elevation  above  the  datum  elevation  of  the  floor, 
and  the  upper-layer/lower-layer  interface;  and  the  pressure  at  the  floor  above  the  specified  datum 
pressure.  Specified  in  an  outside  space  are:  a reference  elevation  above  the  datum  elevation,  and  the 
pressure  at  this  reference  elevation  above  the  specified  datum  pressure. 

When  the  upper  gas  is  less  dense  than  the  lower  gas,  i.e.,  the  fluid  configuration  is  stable,  the  flow 
through  the  vent  Is  determined  by  a traditional  orifice-type  flow  model.  Then,  flow  is  determined  by  the 
cross-vent  pressure  difference  without  any  regard  for  buoyancy  effects  (see,  e.g.,  [2]  and  [3])  and  the 
present  algorithm/subroutine  is  Identical  to  that  of  VENTCF  [4].  When  the  configuration  is  unstable  and 
the  upper  gas  is  more  dense  than  the  lower  gas,  the  effects  of  combined  pressure  and  buoyancy  forces 
can  be  significant.  For  example  if  the  cross-vent  pressure  is  relatively  small,  the  unstable  density 
configuration  leads  to  an  exchange-type  of  flow,  with  gas  in  the  lower  space  rising  into  the  upper  space 
and  gas  from  the  upper  space  dropping  into  the  lower  space,  where  the  flow  rate  from  the  high-  to  the 
low-pressure  side  of  the  vent  is  the  larger  of  the  two.  Also,  even  when  the  cross-vent  pressure  difference 
large  enough  to  produce  uni-directional  flow,  the  effect  of  buoyancy  can  be  great  enough  to  reduce 
significantly  the  flow  rate  from  what  it  would  be  in  the  absence  of  a cross-vent  density  difference.  In  the 
present  algorithm,  for  the  case  of  unstable  configurations  the  calculation  of  the  flow  between  the  two 
spaces  is  based  on  the  theory  developed  and  presented  in  [5]. 

For  unstable  density  configurations,  the  model  of  [5]  and,  therefore,  the  VENTCF2  and  VENTCF2A 
algorithms/subroutines  are  for  flow  through  a circular,  shallow  (i.e.,  small  ratio  of  depth  to  diameter), 
horizontal  vent,  it  is  also  expected  that  the  model  will  give  reasonable  estimates  of  flow  even  for  non- 
circular  vents,  provided  the  aspect  ratio  (maximum-to-minimum  span)  of  a vent  shape  of  interest  is  not 
too  much  different  than  1.  Indeed,  one  of  the  example  calculations  of  [5],  which  includes  comparisons 
with  some  relevant  experimental  data,  provides  limited  support  for  the  applicability  of  the  model  in  the 
case  of  square  vents.  However,  USE  OF  VENTCF2  AND  VENTCF2A  IN  HIGH-ASPECT-RATIO  AND/OR 
MODERATE-TO-LARGE-DEPTH  VENT  SCENARIOS  IS  NOT  VALID  when  cross-vent  pressure  differences 
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are  small-to-moderate  compared  to  App,^ooQ.  Results  beyond  those  developed  in  [5]  are  required  before 
the  present  work  can  be  used  for  the  latter  types  of  vent  shape. 

The  geometry  and  the  conditions  local  to  the  vent  which  determine  the  characteristics  of  the  vent  flow  are 
depicted  in  Figure  2.  These  Include:  the  densities,  and  p^,  and  the  hydrostatic  pressures,  p^  and  pg, 
at  the  elevation,  but  away  from  the  immediate  vicinity  of  the  vent  in  the  upper  and  lower  spaces, 
respectively,  and  the  area,  A^,  and  shape  of  the  vent.  Regarding  the  shape,  at  the  present  time  results 
for  horizontal  vent  flows  are  only  available  for  circular  or  square  vents.  Other  properties  local  to  the  vent 
and  indicated  In  Figure  2 are  T^  and  Tg,  the  absolute  temperatures,  Cog  ^ and  Cog  g,  the  concentrations  of 
oxygen,  and  c^^  and  c^g,  K = 2 to  Nppojj,  the  concentrations  of  a product  of  combustion  K. 

To  simulate  fire  scenarios  at  times  when  vent  flows  are  driven  by  arbitrarily  high  cross-vent  pressure 
differences  (I.e.,  when  compressibility  effects  begin  to  be  significant),  whether  for  stable  or  unstable  cross- 
vent density  configurations,  both  VENTCF2  and  VENTCF2A  implement  the  Ideas  introduced  in  reference 
[3]  and  Implemented  previously  in  the  algorithm/subroutines  VENTHP  [6]  and  VENTCF  [4].  Such  high 
cross-vent  pressure  differences  could  occur,  for  example,  in  fire  scenarios  Involving  flows  through  cracks 
In  otherwise  hermetically-sealed  fire  compartments. 

The  feature  which  defines  the  diffence  between  VENTCF2  and  VENTCF2A  is  in  the  treatment  of  the  source 
of  the  vent  flow  when  the  adjacent-vent  layer  in  the  source  room  is  relatively  thin.  When  the  thickness 
of  the  vent-adjacent  layer,  6y,  is  large  enough,  it  is  assumed  in  VENTCF2A  that  the  flow  driven  through 
the  vent  is  extracted  entirely  from  this  layer.  Unless  5y  = 0,  this  assumption  is  always  made  in  the 
VENTCL2  model/algorithm,  i.e.,  independent  of  the  size  of  6y.  However,  in  VENTCL2A,  when  6y  is  smaller 
than  some  specified,  non-zero,  characteristic  thickness,  6,  the  model/algorithm  assumes  that  only  the 
fraction  6y/5  of  the  net  volumetric  rate  of  flow  driven  through  the  vent  from  the  room  is  extracted  from  the 
vent-adjacent  layer,  the  remaining  fraction,  (1  - 6y/6),  being  extracted  from  the  far  layer. 

It  is  reasonable  to  take  the  characteristic  thickness,  6,  as  being  of  the  order  of  the  characteristic  span  of 
the  vent  opening,  taken  as  to  equivalent  vent  diameter,  D = (4Ay/Tr)^^,  but  never  greater  than  a 
characteristic  distance  related  to  the  height  of  the  room,  H,  from  which  the  flow  is  extracted.  Here  6 is 
taken  to  be  6 = min(D/2,  H/2). 


OUTPUT 

^02,1  [^02.2! 

Concentration  of  oxygen  of  the  flow  from  the  top  [bottom]  space  as  it  enters  the  vent,  [(kg  of 
oxygen)/(kg  of  layer)] 


^K,1  [^K.2] 

Concentration  of  product  K of  the  flow  from  the  top  [bottom]  space  as  it  enters  the  vent,  [(unit 
of  product  K)/(kg  of  layer)] 

^NT.02.p  I = 1 or  2 

Concentration  of  oxygen  in  the  component  of  the  vent  flow  entering  space  I,  provided  such  vent 
flow  component  is  non-zero,  I.e.,  provided  is  non-zero,  [(kg  of  oxygen)/(kg  of  vent  flow)] 


^NT.K.iI  I = 1 or  2;  K = 2 to  NPROD 
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Concentration  of  product  of  combustion  K in  the  component  of  the  vent  flow  entering  space  I, 
provided  such  vent  flow  component  is  nonzero,  i.e.,  provided  is  non-zero,  [(unit  of 

product)/(kg  of  vent  flow)] 

^U.l  I = 1 Of*  2 

If  space  I is  an  inside  room:  Rate  at  which  mass  is  added  to  the  upper  [lower]  layer  of  room  I 
due  to  the  vent  flow  component  which  enters  the  other  space.  Note  that  this  will  always  be 
negative,  since  the  layer  which  supplies  the  material  flowing  to  the  other  space  will  always  have 
mass  extracted  from  It.  [W] 

If  space  I is  an  "outside  room":  | is  the  rate  at  which  mass  is  added  to  space  I due  to  the  vent 
flow  component  which  enters  the  other  space  (i.e.,  the  inside  room).  Note  that  this  will  always 
be  negative,  since  an  outside  space  which  supplies  material  flowing  through  the  vent  to  an 
adjacent  space  will  always  have  mass  extracted  from  it.  is  identical  to  My ,.  [W] 

^vENT,i>  I = 1 to  2 


Mass  flow  rate  of  the  vent  flow  component  entering  space  I.  [kg/s] 

^02.Li  l^02.u.il  i I = or  2 

If  space  I is  an  inside  room:  Rate  at  which  oxygen  is  added  to  the  upper  [lower]  layer  of  room 
I due  to  the  vent  flow  component  which  enters  the  other  space.  Note  that  this  will  always  be 
negative,  since  the  layer  which  supplies  the  material  flowing  through  the  vent  to  the  other  space 
will  always  have  its  convected  oxygen  extracted  from  it.  [(kg  of  02)/s] 

If  space  I is  an  "outside  room":  ^ which  oxygen  is  added  to  space  I due  to  the 

vent  flow  component  which  enters  the  other  space  (i.e.,  the  inside  room).  Note  that  this  will 
always  be  negative,  since  an  outside  space  which  supplies  material  flowing  through  the  vent  to 
an  adjacent  space  will  always  have  its  convected  oxygen  extracted  from  it.  P02.U1  identical  to 

^02,U,l*  02)/®] 

^02.VENT.I>  I = 1 to  2 

Mass  flow  rate  of  oxygen  of  the  vent  flow  component  entering  space  I.  [(kg  of  oxygen)/s] 

^K.L,i  I^K.u.il'  K = 2 to  NppQp;  I = 1 or  2 

if  space  I is  an  inside  room:  Rate  at  which  product  of  combustion  K is  added  to  the  upper  [lower] 
layer  of  room  I due  to  the  vent  flow  component  which  enters  the  other  space.  Note  that  this  will 
always  be  negative  since  the  layer  which  supplies  the  material  flowing  through  the  vent  to  the 
other  space  will  always  have  its  convected  product  of  combustion  K extracted  from  it.  [(unit  of 
product  K)/s] 

If  space  I Is  an  "outside  room":  R,^y , is  the  rate  at  which  product  of  combustion  K is  added  to 
space  I due  to  the  vent  flow  component  which  enters  the  other  space  (i.e.,  the  inside  room).  Note 
that  this  will  always  be  negative,  since  an  outside  space  which  supplies  material  flowing  through 
the  vent  to  an  adjacent  space  will  always  have  the  convected  product  of  combustion  K extracted 
from  it.  Pk  1 1 is  Identical  to  P^ui-  [(unit  of  product  K)/s] 

Pk,vent.,:  I = 1 to  2;  K = 2 to  Np„oD 
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Flow  rate  of  product  of  combustion  K in  the  vent  flow  component  entering  space  I.  [(unit  of 
product  K)/s] 


p,  [P2I 


Absolute  hydrostatic  pressure  in  the  top  [bottom]  space  at  the  elevation  of  the  vent.  [Pa  = 
kg/(m-s^] 

^u,i  I = 1 or  2 

If  space  I is  an  inside  room:  Rate  at  which  enthalpy  added  to  the  upper  [lower]  layer  of  room  I 
due  to  the  vent  flow  component  which  enters  the  other  space.  Note  that  this  will  always  be 
negative,  since  the  layer  which  supplies  the  material  which  flows  to  the  other  space  will  always 
have  its  convected  enthalpy  extracted  from  it.  The  enthalpy  is  based  on  the  absolute  temperature 
of  the  flow,  Tvent.i-  [W] 

If  space  I is  an  "outside  room":  Oy , is  the  rate  at  which  enthalpy  is  added  to  space  I due  to  the 
vent  flow  component  which  enters  the  other  space  (i.e.,  the  inside  room).  Note  that  this  will 
always  be  negative,  since  an  outside  space  which  supplies  material  which  flows  through  the  vent 
to  an  adjacent  space  will  always  have  its  convected  enthalpy  extracted  from  It.  The  enthalpy  is 
based  on  the  absolute  temperature  of  the  flow,  is  Identical  to  Oy ,.  [W] 

^vENT.i*  I = 1 to  2 

Total  enthalpy  flow  rate  in  the  vent  flow  component  entering  space  I.  This  is  based  on  the 
absolute  temperature  of  the  flow,  Ty^^,.  [W] 


Ti[TJ 


Absolute  temperature  of  the  flow  from  the  top  [bottom]  space  as  it  enters  the  vent.  [K] 
**^vENT.i*  I = 1 to  2 


Absolute  temperature  of  the  vent  flow  component  entering  space  I,  provided  such  vent  flow 
component  is  non-zero.  [K] 


Ap 


P2  - p.,,  i.e.,  difference  between  the  pressure  in  the  bottom  space,  at  the  elevation  of  the  vent,  and 
the  pressure  in  the  top  space,  at  the  elevation  of  the  vent.  [Pa  = kg/(m  • s^] 


Density  of  the  flow  from  the  top  [bottom]  space  as  it  enters  the  vent,  [kg/mg] 

^vENT.i*  I = to  2 

Density  of  the  vent  flow  component  entering  space  I,  provided  such  vent  flow  component  is 
nonzero,  [kg/m^] 


INPUT 
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Av 


Area  of  the  vent  [m^]. 

^LK.1  K = 2 to  NppQQ,  I = 1 or  2 

If  space  I Is  an  Inside  room:  Concentration  of  product  of  combustion  K in  lower  [upper]  layer  of 
room  I If  the  volume  of  the  lower  [upper]  layer  is  non-zero  (if  the  lower  [upper]  layer  volume  is 
zero  then  the  c^^,  [c^  ^,]  value  Is  not  used  In  the  calculation),  [(unit  of  product  K)/(kg  of  layer)] 

If  space  I Is  an  ■outside  room':  c^^.  Is  the  uniform  concentration  throughout  the  space  of  product 
of  combustion  K;  c^  k.i  is  specified  as  being  identical  to  c^^i.  [(unit  of  product  K)/(kg  of  local 
atmosphere)] 

^L,02.l  I^U,02.|1'  I = "i  or  2 

If  space  I is  an  inside  room:  Concentration  of  oxygen  in  lower  [upper]  layer  of  room  I if  the 
volume  of  the  lower  [upper]  layer  is  non-zero  (if  the  lower  [upper]  layer  volume  is  zero  then  the 
c^02.i  ICu.o2,il  i®  calculation),  [(kg  of  oxygen)/(kg  of  layer)] 

If  space  I is  an  'outside  room':  c^og.!  i®  uniform  concentration  throughout  the  space  of 
oxygen;  c^  q2.i  is  specified  as  being  Identical  to  c^02.i-  [(^9  of  oxygen)/(kg  of  local  atmosphere)] 

Cp 

specific  heat  at  constant  pressure  of  the  vent  flow.  [W  • s/(kg  • K)]  (suggest  1 0^  W • s/(kg  • K)  for  air 
as  default) 

i^PMAX 

Maximum  allowed  number  of  products  of  combustion. 

^PROD 

Number  of  products  of  combustion,  including  oxygen,  being  tracked  in  the  simulation. 


Datum  absolute  pressure.  [Pa  = kg/(m  • s^] 

"^Li  I = 1 or  2 

If  space  I is  an  inside  room:  Absolute  temperature  of  the  lower  [upper]  layer  in  room  I if  the 
volume  of  the  lower  [upper]  layer  is  non-zero  (if  the  lower  [upper]  layer  volume  is  zero  then  the 
\i  rTu.il  ''^ioe  is  not  used  in  the  calculation).  [K] 

If  space  I is  an  'outside  room':  T^,  is  the  uniform  absolute  temperature  there,  taken  to  be  the 
temperature  at  the  reference  elevation,  ypEpii  , is  specified  as  being  identical  to  T^,.  [K] 


ycEiL.1  [yREF.il 
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If  space  I is  an  inside  room:  Elevation  of  the  ceiling  [floor]  of  room  I above  the  datum  elevation, 
[m] 


If  space  I is  an  'outside  room':  y^EiL.!  and  yp^p,  are  both  identical  and  equal  to  the  reference 
elevation  of  space  I above  the  datum  elevation,  i.e.,  the  specification  must  satisfy  ycpm,  = yppp,. 
The  latter  Identity,  which  will  never  be  satisfied  for  an  inside  room,  is  a characteristic  of  the  input 
data  used  to  distinguish  an  inside  room  from  an  'outside  room.'  [m] 

VLAYEai’  I = "I  01'  2 

If  space  I is  an  inside  room:  Elevation  of  the  upper/lower  layer  interface  in  room  I above  the 
datum  elevation,  [m] 

If  space  I Is  an  'outside  room':  yLAYEai  specified  as  being  identical  to  yppp,.  [m] 


VvENT 


Elevation  of  the  vent  above  the  datum  elevation.  Note  that  y^^,^  must  be  identical  to  either  ycE\u 
or  yppp,  for  each  of  the  one  or  two  inside  rooms  involved  In  the  calculation,  [m] 

^Preri*  I = 1 or  2 

Pressure  at  the  reference  elevation,  yppp,,  in  space  I above  the  datum  absolute  pressure,  Pq^j. 
If  space  I is  an  inside  room,  then  6pppp,  and  yppp,  must  correspond  to  the  pressure  and  elevation, 
respectively,  within  the  room  and  at  the  floor.  [Pa  = kg/(m  • s^] 


e 


p 


Error  tolerance  for  6pppp|.  If  PeRRoai  defined  as  the  uncertainty  In  5pppp,,  I = 1 or  2,  then 
Perrom  satisfies 

IPerror.iI  SPo  I^PRERlI^ 

where  Po  = 1 .OPa  The  first  term  is  based  on  an  absolute  error  tolerance  and  dominates  the 
above  error  bound  when  |6pR^p,|  is  less  than  t.OPa.  The  second  term  is  a relative  error 
tolerance  and  dominates  when  |6pp^p,|  is  greater  than  1.0Pa.  should  be  chosen  to  be 
consistent  with  the  tolerance  specified  for  the  computation  of  6pR£p,  terms  in  the  overall 
compartment  fire  model  computer  code  which  uses  this  algorithm. 

^ui  I = 1 or  2 

If  space  I is  an  inside  room:  Density  of  the  lower  [upper]  layer  in  room  I if  the  volume  of  the  lower 
[upper]  layer  is  non-zero  (if  the  lower  [upper]  layer  volume  is  zero  then  the  [py  J value  is  not 
used  in  the  calculation),  [kg/m^] 

If  space  I is  an  'outside  room':  p^  is  the  uniform  density  there;  py , is  specified  as  being  identical 
to  pl,,.  [kg/m^] 


CALCULATION 
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Set  all  output  to  zero.  If  = 0 then  calculation  is  complete,  if  not  then: 


1.  Calculate  p,  for  I = 1 and  2 , and  then  Ap  (the  p,  calculation  follows  the  DELP  algo- 
rithm/subroutine of  reference  [7]): 


Pi  “ ‘^Pref.i  + ^Pi  Pdat 

^P  P2  ■ Pi  ~ (^PrEF,2  ■ ^PrEF.i)  (^P2  ■ ^Pl) 


(1) 


(2) 


where 


VrEF.I  - VvENT  - yiAYER,!  ^^eni  6P|  - - Pmg(yvENT  ■ Vref.i) 
else:  6p|  = - PL,ig(yLAYER.i " Vref.i)  ’ ^u.iQCYvent  " Vlayer,!) 


(3) 


and  where  g,  the  acceleration  of  gravity,  is  9.8  m/s^. 


2.  Calculate  D,  6,,  6y,,  p^,  T,,  c^gj.  and  the  c,^,.  K = 2 to  NpRop,  for  I 


PROD' 


for  I = 1 and  2,  and  then  Ap: 


D = (4Av/71)’'^ 


(4) 


VcEiui  ^ Vref.i  (space  I is  an  Inside  space  and) 

6|  = min[(ycEiL,|  - yREF,i)/2i  D/2];  6y,  = |yLAYER,i " YventI 
else  (Yceiui  ~ Vref.i'  apace  I is  an  outside  space,  and) 

6,  = D/2;  6y,  = 0 

If  6y,  > 6,  then  (vent  flow  from  space  I is  extracted  from  the  vent-adjacent  layer): 
If  I = 1 (i.e.,  the  upper  room)  then 


else  (I  = 2,  i.e.,  the  lower  room  and) 


else  (6y,  < 6,  and  the  vent  flow  from  space  I is  extracted  from  both  layers) 

If  I = 1 (I.e.,  the  upper  room)  then 

P-\  ~ "^1  ~ 

^02.1  “ I^U,02.1^U,l(^  ■ ^Yi/^l)  ■f’  ^L,02.1^L,1^y/^ll/^1' 

^K.1  ~ " ^yi/^i)  ^ ~ 2,  NppQQ 

else  (I  = 2,  i.e.,  the  lower  room  and) 
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P2  ~ Pl,2^^  " Pu,2^y2^^2’  "^2  ~ '^L,2^L,2^^2’ 

^02,2  ~ I^L,02,2^L,2(^  ~ ^)l2!^^  ^\i,02,2P\i,2^y 2!^:^! P2' 

^K,2  ~ I^L,K2^L2(^  " ^U.K2P\J,2^y2^^2y P2*  K = 2,  NppQQ 


= Pl  - P2  (5) 

Define  Vsjhigh  volume  rate  of  flow  through  the  vent,  from  the  high-  to  the  low-pressure 
space,  that  is  predicted  with  a "standard,”  unidirectional-flow-type  calculation  (i.e.,  without  regard 
to  the  effect  of  buoyancy,  in  general,  or  the  stability  of  the  cross-vent  density  configuration,  in 
particular),  where  arbitrarily  high  cross-vent  pressures  are  allowed.  Here,  the  calculation  follows 
the  model  of  [3]  as  implemented  In  the  VENTHP  algorithm/subroutine  of  [6].  Calculate  Vst.high* 

If  Ap  = 0 then:  VgjniGH  = ^ CALCULATION  (6) 

If  Ap  > 0 or  Ap  < 0 then: 


a. 


b. 


c. 


Define  and  compute  phigh* 

If  Ap  > 0,  i.e.,  "standard"  flow  from  (lower)  space  2 to  (upper)  space  1,  then: 


PhIGH  "■  P2»  ^ ^p/p2 

If  Ap  < 0,  i.e.,  "standard"  flow  from  (upper)  space  1 to  (lower)  space  2,  then: 

Phigh  “ Pi»  ^ ~ ^P/Pi  (8) 

X = 1 - e (9) 

Compute  C(x),  the  vent  flow  coeffient,  and  w(x): 

C(x)  = 0.85  - 0.25X  = 0.60  + 0.25e  (10) 

' 1 - [3/(47)]e  if  0 < £ < 10  ® 

w(x)  = - (11) 

f(x)/[2e]^^  if  1 > e > 10  ® 

where 


f(x)  = 


^ {[27/(7  - 1)]x^[1  - if  e < 1 - [2/(7  + 1)r'^  ■ 

, {712/(7  + 1)]*^  if  E a 1 - [3/(7  + • ’• 


(12) 


and  where  7,  the  ratio  of  specific  heats  of  the  vent  flow  gas,  is  taken  to  be  that  of  air, 
1.40.  Note  that  for  the  present  horizontal  vent  application,  the  C(x)  of  Eq.  (10)  is  taken 
to  be  consistent  with  the  standard  incompressible  limit  for  flow  through  circular  sharp- 
edged  orifices  in  the  sense  that  C 0.60  as  e -►  0 [8]. 

Define  and  compute  Ap^^: 

~ [epM/0((i.OPa,  |6p„EP,|.  (13) 
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d.  Define  and  compute  a numerical  damping  factor,  and  then  Vsjhigh- 

^NoisE  ~ 1.0  - exp(-|Ap|  ^/APq^  (14) 

^ST.HIGH  ~ ^NOISE^W'^W(^^HIGh)  ^^vI^PI  ^ 

The  term  F^ojg^  of  Eq.  (14)  is  designed  to  damp  out  the  numerical  noise  (error)  in  the 
calculated  value  for  Ap  that  would  otherwise  be  dominant  in  Eq.  (1 5)  when  Ap  is  small 
relative  to  the  maximum  of  1.0  pascal  and  the  calculated  reference  pressures,  6pp^p^, 
5Prep2.  The  term  Ap^ur  of  Eq.  (13)  Is  an  estimate  of  how  small  the  maximum  of  |Ap| 
must  be  to  retain  a few  digits  of  accuracy  in  the  calculation  of  Ap.  When  the  calculated 
value  of  I Ap|  Is  smaller  than  p^jj,  this  value  and,  therefore,  the  value  of  VgjRiQR  in  Eq. 
(15)  will  likely  contain  noise  which  should  be  damped.  F^Qjgp  is  constructed  to  tend 
towards  1 when  | Ap|  is  large  relative  to  Ap^uj  and  tends  towards  0 when  | Ap|  is  small 
relative  to  Ap^uj. 

When  the  cross-vent  density  configuration  is  unstable  {Ap  > 0),  mixed  pressure-  and  buoyancy- 
driven  aspects  of  the  flow  have  to  be  considered  and  the  vent  flow  rates  are  calculated  according 
to  [5].  Define  Vb^iqr,  Vg  Low  reference- [5],  buoyancy-affected,  volume  flow  rates  across 
the  vent  from  the  high-to-low  and  low-to-high  pressure  spaces,  respectively. 

If  Ap  < 0:  Vb  high  = ^B.Low  = skip  to  (the  next)  Item  5 of  the  CALCULATION  (17) 

a.  Calculate  f,  p,  /i(T)  In  m^/s  from  [12],  and  e > 0 : 

T = (T,  + T2)/2;  p = (p,  + P2)I2-.  p(f)  = p[0.04128(10-')T“/(f  + 110.4)); 

£ = Ap/p  >0  (18) 

b.  Calculate  Gr;  if  Pg  > p^,  i.e.,  Ap  > 0,  then  replace  £ by  - e < 0: 

Gr  = 2gD®|£l/[/i(T)/p]^;  if  Ap  > 0 then  £ = - £ (19) 

c.  Calculate  conditions  at  the  limit  of  unidirectional  flow  (i.e.,  the  flooding  condition)  and  the 
relative  pressure,  6p*: 

^HiGH, FLOOD  ~ 0.1754AY(2gD j£|)^^exp(0.5536£); 

^Pflood  = 0.2427(4gApD)(1  + £/2)exp(1.1072£);  (20) 

5p*  = I Ap|/App|_QGQ 

d.  Calculate  a^,  Og  and  then  Vq  low  Vb.high- 

= Fno,seC(x)w(x)/0.1780;  = 1.045  (21) 

If  6p*  > 1 , exp)ect  uni-directional  flow: 

^B.LOW  “ ^B.HIGH  ” ^HIGH.FLOOD'C^  " ^2^  (22) 


If  6p*  < 1 , expect  mixed  flow: 
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^EX,MAX  = 0.055(4/7r)Av(gD|e|)'/2.  ^ . q jqjq.  ^ ^ ^23) 

^B.LOW  ~ ^EXMAxK^  1^3/2) (1  - 5p*)  - (2  + 1X13/2) (1  - 6p*)] 

(24) 

^B.HIGH  ~ ^B.LOW  - [1  + (M  - 1)(1  - 6p*)]  ^}^hIGH,FLOOc/(^  " ^) 

(Note:  of  Eq.  (21)  is  somewhat  different  than  in  [5].  The  modification  here  is 

provide  for  an  analytic  representation  of  vent  flow  which  is  continuous  and  uniformly  valid 
even  as  5p*  is  increased  to  a level  where  compressibility  effects  become  important.  The 
result  of  [5]  does  not  Include  the  effect  of  compressibility.) 

5.  Following  [SJ^  when  Ap  > 0 the  above  will  be  taken  as  the  flow  solution  provided  Gr  > 2(10^). 
For  smaller  Gr,  the  reference-[5]  sojution  begins  to  loose  Its  validity,  and  there  is  no  existing 
model  for  the  flows  the  range  0 < Gr  <_2(10^).  However,  it  Is  clear  that  the  above-calculated 
“standard'  flow  must  be  approached  as  Gr  -►  0. 

Consistent  with  the  above  comments,  when  Gr  is  in  the  range  0 < Gr  < 2(10^),  the  vent  flow  will 
be  estimated  by: 

flow  = (“standard"  flow)  -i-  [(reference-[5]  flow)  - (“standard"  flow)]Gr/[2(10^)] 

a.  Define  and  Vlqw  volume  flow  rates  across  the  vent  from  the  high-to-low  and  low- 
to-high  pressure  spaces,  respectively. 


if  Ap 

- 9-  ^HIGH  ” ^ST.HIGH*  ^LOW 

= 0 

(25) 

If  Ap 

> 0 and  0 < Gr  < 2(10^): 

^HIGH  ” ^ST.HIGH  (^B.HIGH 

■ ^ST. high)® *^^[20  9 )] 

(26) 

^LOW  ” ^b.low®*^^[2(1  0 )] 

If  Ap 

> 0 and  Gr  > 2(10^): 

” ^B.HIGH'  ^LOW  “ ^B.LOW 

(27) 

b.  Define  and  calculate  Vy^Ni.i,  I = 1 and  2,  as  the  vent  volume  flow  rate  entering  space  I. 

If  Ap  > 0:  = ^HIGH>  ^VENT,2  “ ^LOW  (2®) 

(i.e.,  flow  to  the  upper  space  is  the  flow  from  the  high-  to  the  low-pressure  space,  and 
flow  to  the  lower  space  is  the  flow  from  the  low-  to  the  high-pressure  space) 

If  Ap  < 0:  ^VENT.2  “ ^HIGH  (29) 

(i.e.,  flow  to  the  upper  space  is  the  flow  from  the  low-  to  the  high-pressure  space,  and 
flow  to  the  lower  space  is  the  flow  from  the  high-  to  the  low-pressure  space) 

0.  Calculate  the  vent  flow  properties  *^venti>  ^^ento2.I’  ^vent,k,I’  I^  ~ 2 to  HppQQ,  I = 1 and  2. 


^VENT.I  “ P2(Pl/p2)'  ^VENT,2  “ P^iP2^P^)•  '^’vENT.1  “ ^2*  '*^VENT.2  “ ^1*  (®®) 
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^ENT,02.1  “ ^02.2'  ^ENT,02.2  “ ^02.1*  ^ENT,K.1  ” ^K.2'  ^ENT,K.2  “ ^K.1'  K - 2 tO  NppQQ 

Calculate  the  vent  flow  rates  I = 1 and  2: 

^VENTI  “ PvENT.rVENT.I  (3"*) 

Calculate  the  vent  flow  rates  02venti’  ^ kventi*  ” 2 to  NppQQ,  I — 1 and  2^ 

^VENT.I  “ ^VENT.I^p'^VENTI'  ^02.VENT.I  “ ^VENT,I^ENT.02.I’ 

(32) 

^KVENT.I  ” ^VENT.I^ENT.K,!’  K = 2 tO  NppQQ 

Calculate  rates  at  which  flows  are  added  to  layers  of  each  space  as  a result  of  the  vent  flow 
extracted  from  it.  First  consider  space  1 (I  = 1 , J = 2)  and  then  space  2 (I  = 2,  J = 1).  For  either 
case: 

If  6y,  > 6|  then  (vent  flow  from  space  I is  extracted  from  the  vent-adjacent  layer): 

If  I = 1 (I.e.,  the  upper  room)  then 

^L.1  ” ■ ^VENT.2«  ^U.l  “ ^L.1  “ ’ ^VEMT.2'  ^U.1  “ 

^02,1,1  “ ■ ^02.VENT,2'  ^02.U.1  ~ 

^K.L.1  ~ ” ^K.VENT.2«  ^K.U.1  = Oi  K = 2,  NppQ^ 

(33) 

VcEiui  ” yREF.1  (room  1 is  an  outside  room  and) 

^U.l  “ ^U.1  “ ^L.1' 

^02.U,1  ” ^02.U1’  ^K.u.1  “ ^K.L.1'  K = 2,  NppQp 
else  (I  = 2,  i.e.,  the  lower  room  and) 

Mu2  = ’ ^VENT.V  ^L.2  “ ^U.2  ” ‘ ^VENT.1'  ^L,2  “ 

^02.U.2  “ ■ ^02,VENT.1'  ^02.L2  “ 

^K.U.2  “ " ^K.VENT.1'  ^KL.2  = 0.  K = 2,  NppQp 

(34) 

ycEiU2  = yREF.2  (room  2 Is  an  outside  room  and) 

Mu  2 = ^L.2»  ^U.2  “ ^U2' 

^02.U.2  ~ ^02.U2'  *^K.U.2  “ ^K.L,2'  ^ = 2,  NppQp 

else  (5y,  < 6,  and  the  vent  flow  from  space  I is  extracted  from  both  layers) 

If  I = 1 (i.e.,  the  upper  room)  then 

^u,i  “ ■ ■ ^yi/^i)^vEm;2'  ^ui  ” ■ ^L,i(^yi/^i)^vENT,i 
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(35) 


^02,U.1  ~ ^U.1^U.02.1'  ^02,L,1  “ ^L.1^U02.1 
^K.U.1  “ ^U,1^U,K.1’  ^K.L,1  “ ^L,1^L.K,1«  ^ = 2,  Np^QQ 
VcEiLi  Vref.i  (room  1 is  an  outside  room  and) 

^L.1  = ^U.1*  ^ui  “ ^U.1» 

^02.U1  ~ ^02,U.1»  ^K,L.1  “ ^K,U.1'  ^ = 2,  NppQ^ 

else  (I  = 2,  i.e.,  the  lower  room  and) 

^U2  “ " Pl,2^^  " ^U.2  “ " 

^U2  “ ^L.2^p*'^U2*  ^U.2  “ ^U.21^p^U.2' 

^02,1,2  “ ^l,2^U02.2'  ^02.U,2  ~ ^U.2^U.02.2« 

^K.L.2  “ ^L2^L,K.2«  ^K,U,2  “ ^U.2^U.K.2’  K = 2,  Np^Q^;  (36) 

ycEiu2  = yREF.2  (room  2 is  an  outside  room  and) 

My  2 ” ^L,2'  ^U,2  “ ^L,2* 

^02.U.2  ” ^02,L,2’  ^K.U.2  “ ^K.L,2'  K = 2,  Np^Qp 
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SUBROUTINE  VARIABLES 

All  nomenclature  In  the  subroutine  is  Identical  to  the  nomenclature  used  above  except  for: 


Ay 

- 

AVENT  [m^] 

C 

- 

COEF 

Cp 

- 

CP  [W.s/(kg-K)] 

^L.KI’  ^U.K.1 

- 

CONL(K,l),  CONU(K,l),  1 = 1 or  2 [(unit  of  product  K)/(kg 
of  layer)] 

^U02.l'  ^U.02.1 

- 

CONL(1,l),  CONU(1,l),  1 = 1 or  2 [(kg  of  oxygen)/(kg  of 
layer)] 

^NT.K,! 

- 

CVENT(K,I),  1=  1 or  2 [(unit  of  product  K)/(kg  of  vent 
flow) 

^NT.02.1 

- 

CVENT(1,I),  1 = 1 or  2 [(kg  of  oxygen)/(kg  of  vent  flow)] 

Ck.1 

- 

C(K,I),  1=  1 or  2 [(unit  of  product  K)/(kg  of  vent  flow) 

^02,1 

" 

C(1,l),  1 = 1 or  2 [(kg  of  oxygen)/(kg  of  vent  flow)] 

^NOISE 

- 

FNOISE  [dimensionless] 

f 

- 

FF  [dimensionless] 

Gr 

- 

GR  [dimensionless] 

g 

- 

9.8  m/s^ 

M 

- 

XM  [dimensionless] 

My  , 

- 

XML(I),  XMU(I),  1 = 1 or  2 [kg/s] 
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^VENT.I 

XMVENT(I),  1 = 1 or  2 [kg/s] 

m3 

XM3  [dimensionless] 

NpMAX 

NPMAX  [dimensionless] 

^PROD 

NPROD  [dimensionless] 

P P 

^K.UI’  ^K.U,I 

PL(K,I),  PU(K,I),  I = 1 or  2 [(unit  of  product  K)/s] 

^K.VENT.1 

PVENT(K,I),  I = 1 or  2 [(unit  of  product  K)/s] 

P P 

^02.UI’  ^02.U.I 

PL(1,I),  PU(1,I),  I = 1 or  2 [(kg  of  oxygen)/s] 

b 

^02.VEMT.I 

PVENT(1,I),  I = 1 or  2 [(kg  of  oxygen)/s] 

Pdat 

PDATUM  [Pa  = kg/(m  • s^] 

Pi 

P(l),  I = 1 or  2 [Pa  = kg/(m-s^] 

0^1,  Ou , 

QL(I),  QU(I),  I = 1 or  2 [W] 

^vE^^;l 

QVENT(I),  I = 1 or  2 [W] 

T 

TBAR  [K] 

T, 

T(l),  I = 1 or  2 [K] 

"^ui’  "^u.i 

TL(I),  TU(I),  1 = 1 or  2 [K] 

*^VENT,l 

TVENT(I),  1 = 1 or  2 [K] 

VBHIGH  = VB(1)  if  Ap  > 0,  = VB(2)  if  Ap  < 0 [m^/s] 

^B.LOW 

VBLOW  = VB(2)  if  Ap  > 0,  = VB(1)  if  Ap  < 0 [m^/s] 

'^EXMAX 

VEXMAX  [m^/s] 

^LOW 

VHIGH,  VLOW  [m®/s] 

^HIGH.FLOOD 

VHIGHFL  [mVs] 

^ST.HIGH 

V = VST(1)  if  Ap  > 0,  = VST(2)  if  Ap  < 0 [m^/s] 

^VEMT.I 

WENT(I),  1 = 1 or  2 [m^/s] 

W 

W [dimensionless] 

X 

X [dimensionless] 

VcEIUI 

YCEIL(I),  1 = 1 or  2 [m] 
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y LAYER.  I 

yREF.I 

yVENT 

7 


Ap 


^Pflood 


Lp 

<5, 

6p, 

^Pref.i 

^yi 

£ 

Phigh 

p 

P\ 

Pu'  Pu.i 
advent  I 

^2 

PREPARED  BY 

Leonard  Y Cooper 
December  1993 


YLAY(I),  I = 1 or  2 [m] 

YREF(I),  I = 1 or  2 [m] 

YVENT  [m] 

1.40 

DELP  [Pa  = kg/(m.s^] 

DPC1D2  [Pa  = kg/(m-s^l 
DELPFD  [Pa  = kg/(m-s^] 

DELDEN  [kg/m^] 

DEL(I),  I = 1 or  2 [m] 

DPDDPFL  [Pa  = kg/(m-s^] 

DP(I),  I = 1 or  2 [Pa  = kg/(m-s^] 
DPREF(I),  I = 1 or  2 [Pa  = kg/(m  .s^] 
DELY(I),  I = 1 or  2 [m] 

EPS  [dimensionless] 

EPSP  [dimensionless] 

XMEW  [m^/s] 

DENBIGH  [kg/m^] 

DENBAR  [kg/m^] 

DEN(I),  I = 1 or  2 [kg/m®l 
DENL(I),  DENU(I),  I = 1 or  2 [kg/m^] 
DENVNT(I),  I = 1 or  2 [kg/m^] 

SIGMA1  [dimensionless] 

SIGMA2  [dimensionless] 
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(C) 

Figure  1 . The  possible  configurations  of  the  two  spaces  joined  by  a horizontal  ceiling/floor  vent 

with  space  1 above  space  2:  a)  two  inside  rooms;  b)  an  outside  space  above  an 
inside  room;  c)  an  Inside  room  over  an  outside  space. 
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Figure  2.  The  geometry  and  conditions  iocai  to  a horizontai  ceiiing/fioor  vent  which  determine 
the  characteristics  of  the  vent  fiow. 
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oooooooooooooooooooooooooooooooooooooooooooooooo 


SUBROUTINE  VENTCF2A 


SUBROUTINE  VENTCF2A( 

I AVENT,CONUCONU,CP,NPMAX,NPROD,PDATUM,TL,TU, 

I YCEiL,YREF,YLAYYVENT.DPREF,EPSP,DENL,DENU. 

O C,CVENT^UXMU^MVENT,PL,PU,PVENT,P,QL,QU,QVENT,T,TVENT, 

O FRACM.FRACQ.FRACP, 

O DELP.DEN.DENVNT) 

INCLUDE  PRECIS.INC 
C^BEQ 

C***  VENTCF2A  - CALCULATION  OF  THE  FLOW  OF  MASS,  ENTHALPY,  OXYGEN 
C AND  OTHER  PRODUCTS  OF  COMBUSTION  THROUGH  A HORIZONTAL 

C VENT  JOINING  AN  UPPER  SPACE  1 TO  A LOWER  SPACE  2.  THE 

C SUBROUTINE  USES  INPUT  DATA  DESCRIBING  THE  TWO-LAYER 

C ENVIRONMENT  OF  INSIDE  ROOMS  AND  THE  UNIFORM 

C ENVIRONMENT  IN  OUTSIDE  SPACES.  THIS  IS  A VARIATION 

C OF  VENTCF2  WHICH  SMOOTHS  THE  RATE  OF  EXTRACTION  FROM 

C VENT-ADJACENT  LAYERS  WHEN  SUCH  LAYERS  ARE  RELATIVELY 

C THIN. 

C 

C*^  SUBROUTINE  ARGUMENTS 


INPUT 


AVENT 

CONL(K,0 


CONL(K.O 


CP 

NPMAX 

NPROD 

PDATUM 

TL(0 

TU(0 

YCEIL(I) 

YREF(0 

YLAY(0 

YVENT 

DPREF(I) 

EPSP 

DENLO) 

DENU(I) 


OUTPUT 


AREA  OF  THE  VENT  [M*«2] 

CONCENTRATION  OF  PRODUCT  K IN  LOWER  LAYER  OF  AN  INSIDE 
ROOM  I OR  UNIFORM  CONCENTRATION  OF  PRODUCT  K IN  AN 
OUTSIDE  SPACE  I [(UNIT  OF  PRODUCT)/(KO  LAYER)] 
CONCENTRATION  OF  PRODUCT  K IN  UPPER  LAYER  OF  AN  INSIDE 
ROOM  I OR  UNIFORM  CONCENTRATION  OF  PRODUCT  K IN  AN 
OUTSIDE  SPACE  I [(UNIT  OF  PRODUCT)/(KO  LAYER)] 

SPEaFIC  HEAT  [W*S/(KO*K)] 

MAXIMUM  ALLOWED  NUMBER  OF  PRODUCTS 
NUMBER  OF  PRODUCTS  IN  CURRENT  SCENARIO 
DATUM  ABSOLUTE  PRESSURE  [PA  « KO/(M*S**2)] 

TEMPERATURE  OF  LOWER  LAYER  OF  AN  INSIDE  ROOM  I OR 
TEMPERATURE  OF  AN  OUTSIDE  SPACE  I [K] 

TEMPERATURE  OF  UPPER  LAYER  OF  AN  INSIDE  ROOM  I OR 
TEMPERATURE  OF  AN  OUTSIDE  SPACE  I [K] 

HEIGHT  OF  CEIUNG  ABOVE  DATUM  ELEVATION  FOR  AN  INSIDE 
ROOM  I OR  YREF(0  FOR  AN  OUTSIDE  SPACE  I [M] 

HEIGHT  OF  REFERENCE  ELEVATION  FOR  SPACE  I ABOVE  DATUM 
ELEVATION  [M] 

HEIGHT  OF  LAYER  ABOVE  DATUM  ELEVATION  FOR  AN  INSIDE 
ROOM  I OR  YREF(I)  FOR  AN  OUTSIDE  SPACE  I [M] 

HEIGHT  OF  VENT  ABOVE  DATUM  ELEVATION  [M] 

PRESSURE  IN  SPACE  I AT  ITS  REFERENCE  ELEVATION  ABOVE 
DATUM  ABSOLUTE  PRESSURE  [PA  = KG/(M*S*^)] 

ERROR  TOLERANCE  FOR  DPREF  [DIMENSIONLESS] 

DENSITY  OF  LOWER  LAYER  OF  AN  INSIDE  ROOM  I OR  DENSITY 
OF  AN  OUTSIDE  SPACE  I [KO/M**3] 

DENSITY  OF  UPPER  LAYER  OF  AN  INSIDE  ROOM  I OR  DENSITY 
OF  AN  OUTSIDE  SPACE  I [KO/M**3] 


C(K,D 

CVENT(I) 

FRACM(0 

FRACQO) 

FRACP(K,I) 

XML(0 


CONCENTRATION  OF  EACH  PRODUCT  IMMEDIATELY  ABOVE  ON 
SPACE  I - 1)  AND  BELOW  ON  SPACE  I - 2)  THE  VENT 
[(UNIT  OF  PRODUCT)/(KO  LAYER)] 

CONCENTRATION  OF  EACH  PRODUCT  IN  THE  VENT  FLOW 
COMPONENT  ENTERING  SPACE  I [(UNIT  OF  PRODUCT)/(KO  OF 
VENT  FLOW)] 

FRACTION  OF  XMVENTO)  EXTRACTED  FROM  NEAR-VENT  LAYER 
OF  SUPPLY  ROOM 

FRACTION  OF  QVENT(I)  EXTRACTED  FROM  NEAR-VENT  LAYER 
OF  SUPPLY  ROOM 

FRACTION  OF  PVENT(K,0  EXTRACTED  FROM  NEAR-VENT  LAYER 
OF  SUPPLY  ROOM 

RATE  AT  WHICH  MASS  IS  ADDED  TO  THE  LOWER  LAYER  OF  AN 
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XMU(0 

INSIDE  ROOM  1 - 1 0 = 2)  OR  TO  AN  OUTSIDE  SPACE  1 « 
10^2)  DUE  TO  THE  VENT  FLOW  COMPONENT  ENTERING 
SPACE  = 2 0 1)  [KO/S] 

RATE  AT  WHICH  MASS  IS  ADDED  TO  THE  UPPER  LAYER  OF  AN 

XMVENT(I) 

INSIDE  ROOM  1 « 1 0 - 2)  OR  TO  AN  OUTSIDE  SPACE  1 « 1 

0 = 2)  DUE  TO  THE  VENT  FLOW  COMPONENT  ENTERING  SPACE 

1 B 2 0 = 1)  [KG/S] 

MASS  FLOW  RATE  IN  THE  VENT  FLOW  COMPONENT  ENTERING 

PL(K,l) 

SPACE  I [KG/S] 

RATE  AT  WHICH  PRODUCT  K IS  ADDED  TO  THE  LOWER  LAYER  0 

PU(K,0 

AN  INSIDE  ROOM  I = 1 (I  = 2)  OR  TO  AN  OUTSIDE  SPACE  1 
« 1 (1  > 2)  DUE  TO  THE  VENT  FLOW  CONPONENT  ENTERING 
SPACE  1 B 2 (1  « 1)  [(UNIT  OF  PRODUCT)/S] 

RATE  AT  WHICH  PRODUCT  K IS  ADDED  TO  THE  UPPER  LAYER  OI 

PVENT(K,I) 

AN  INSIDE  ROOM  1 = 1 (I  2)  OR  TO  AN  OUTSIDE  SPACE  1 
« 1 (1  « 2)  DUE  TO  THE  VENT  FLOW  COMPONENT  ENTERING 
SPACE  1 - 2 (1  - 1)  [(UNIT  OF  PRODUCT)/S] 

FLOW  RATE  OF  PRODUCT  K IN  THE  VENT  FLOW  COMPONENT 

P(l) 

ENTERING  SPACE  I [(UNIT  OF  PRODUCT)/S] 

ABSOLUTE  PRESSURE  IMMEDIATELY  ABOVE  (IN  SPACE  I > 1) 

QL(I) 

. 

AND  BELOW  (IN  SPACE  1 « 2)  THE  VENT  [PA  > KQ/(M*S*^)] 
RATE  AT  WHICH  ENTHALPY  IS  ADDED  TO  THE  LOWER 

QUO) 

LAYER  OF  AN  INSIDE  ROOM  I = 1 (1  = 2)  OR  TO  AN  OUTSIDE 
SPACE  1 c 1 (1  s 2)  DUE  TO  THE  VENT  FLOW  COMPONENT 
ENTERING  SPACE  1 « 2 (1  [W] 

RATE  AT  WHICH  ENTHALPY  IS  ADDED  TO  THE  UPPER  LAYER  OF 

OVENTO) 

AN  INSIDE  ROOM  1 e 1 (|  . 2)  OR  TO  AN  OUTSIDE  SPACE  1 

B 1 (1  2)  DUE  TO  THE  VENT  FLOW  COMPONENT  ENTERING 

SPACE  1 B 2 (1  b:  1)  [W] 

FLOW  RATE  OF  ENTHALPY  IN  THE  VENT  FLOW  COMPONENT 

TO) 

ENTERING  SPACE  1 [W] 

ABSOLUTE  TEMPERATURE  IMMEDIATELY  ABOVE  ON  SPACE  1 « 

TVENTO) 

1)  AND  BELOW  ON  SPACE  1 « 2)  THE  VENT  [K] 

ABSOLUTE  TEMPERATURE  OF  THE  VENT  FLOW  COMPONENT 

DELP 

ENTERING  SPACE  1 [K] 

CROSS-VENT  PRESSURE  DIFFERENCE.  P(2)  - P(1),  PA  « 

DENO) 

KQ/(M*S**2)] 

DENSITY  IMMEDIATELY  ABOVE  ON  SPACE  1 = 1)  AND  BELOW 

DENVNTO) 

ON  SPACE  1 = 2)  THE  VENT  [KO/M**3] 

DENSITY  OF  THE  VENT  FLOW  COMPONENT  ENTERING  SPACE  1 

[KG/M**3] 

C****  NOTE  THAT  NPMAX2  SHOULD  BE  BlOQER  THAT  NPIMAX 
PARAMETER  (NPMAX2s10) 

DIMENSION  CONL(NPMAX2^).  CONU(NPMAX2^).  TL(2).  TU(2) 
DIMENSION  YCEIL(2),  YREF(2),  YLAY(2),  DPREF(2) 

DIMENSION  DENL(2).  DENU(2) 

DIMENSION  C(NPMAX2^),  CVENT(NPMAX2^).  XML(2),  XMU(2) 
DIMENSION  XMVENT(2) 

DIMENSION  FRACM(2).FRACQ(2),FRACP(NPMAX2^) 

DIMENSION  PL(NPMAX2^),  PU(NPMAX2^),  PVENT(NPMAX2^).  P(2) 
DIMENSION  OL(2),QU(2) 

DIMENSION  QVENT(2),  T(2),  TVENT(2),  DEN(2),  DENVNT(2) 
DIMENSION  DP(2),VB(2),VST(2).  VVENT(2) 

DIMENSION  DEL(2),DELY(2) 

PARAMETER  (QAMs1.4000) 

DATA  IFIRST/0/ 

SAVE  IRRST,QAMCUT,QAMMAX 
PARAMETER  (Q«9.80D0) 

PARAMETER  (PI>3.1415926S4D0) 

PARAMETER  (SIOMA2-1.045D0) 

PARAMETER  (XM3s-0.7070DO) 

OR-O.DO 

APOWER>1.DO 

VHiaHFL-O.DO 

DELPFL-O.DO 
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DO  5 1=1,2 
VST(l)=0.D0 
VB(l)-0.D0 
QVENT(l)=0.D0 
XMVENT(I)=O.DO 
XML(l)=0.D0 
XMU(0«O.DO 
QL(l)«0.D0 
QU(1)>0.D0 
DO  2 NP=1.NPROD 
PVENT(NRI)=O.DO 
PL(NP,l)=0.D0 
PU(NRI)=0.D0 
2 CONTINUE 

5 CONTINUE 

C 

C*^  THE  FOLLOWING  CODE  SEGMENT  COMPUTES  CONSTANTS  REQUIRED  BY  VENTCF2. 
C***  IT  IS  EXECUTED  THE  FIRST  TIME  VENTCF  IS  CALLED. 

C 

iF(IFIRSTEQ.O)THEN 
IFIRST  = 1 

GAMCUT  = (2.0D0/(GAM+1.0D0))**(OAM/(OAM-1.0D0)) 
ZZZ=OAM*((2.0D0/(GAM-l-1.0D0))**((OAM+1.0D0)/(OAM-1.0D0))) 

GAMMAX  = DSQRT(ZZZ) 

ENDIF 

C 

C***  IF  AVENT  = 0 THEN  CALCULATION  IS  OVER 
C 

IF(AVENT.EQ.O.DO)OOTO  200 
C 

C***  1.  CALCULATE  THE  PC)  AND  DELP 
C 

DO  10  I > 1.2 

IF(YREF(l).LE.YVENTJVND.YVENTLE.YLAY(0)THEN 

C 

C***  THE  VENT  IS  AT  OR  BELOW  THE  REFERENCE  ELEVATION  IN  SPACE  I.  IF 

C***  SPACE  I IS  AN  INSIDE  ROOM  THEN  BOTH  THE  VENT  AND  THE  LAYER 

C***  INTERFACE  ARE  AT  THE  FLOOR  ELEVATION. 

C 

DPO  = - 0*DENL(l)*(YVENT  - YREFW) 

ELSE 

C 

C***  THE  VENT  IS  ABOVE  THE  REFERENCE  ELEVATION  IN  SPACE  I. 

C***  IF  SPACE  I IS  AN  INSIDE  ROOM  THEN  THE  VENT  IS  AT  THE 
C***  CEIUNG. 

C 

DP(I)  - • O*0ENL(D*(YLAYO-YREF(l)) 

$ - 0*DENU(l)*(YVEhrPYLAY(l)) 

ENDIF 

PO  - DPREFO  •t'  DPO  + PDATUM 
10  CONTINUE 

C 

C***  DELP  IS  PRESSURE  IMMEDIATELY  BELOW  THE  VENT  LESS  PRESSURE 

C***  IMMEDIATELY  ABOVE  THE  VENT. 

C 

DELP  - pPREF(2)-DPREF(1))-»-pP(2)-DP(1)) 

C 

C***  2.  CALCULATE  D,  DELO.  DELYO  AND  THE  EFFECTIVE  PROPERTIES  OF 

C***  FLOW  INTO  THE  VENT.  IF  AN  ADJACENT^VENT  LAYER  IS  TOO  THIN,” 

C***  THESE  ARE  BASED  ON  MIXING  OF  GAS  EXTRACTED  FROM  BOTH  THE 

C***  ADJACENT  AND  THE  FAR  LAYER.  CALCULATE  DELDEN. 

C 

D«DSQRT(4.DO*AVENT/P0 
DO  31  I - 1,  2 

DENO  - 0.0D0 

TO  - 0.000 

DO  22  K - 1,  NPROD 

C(ig)  - 0.0D0 

22  CONTINUE 
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c 

c*** 

c 


c 

c 


24 


26 


# 


# 

# 


28 


# 


# 

# 


30 


31 

C 

c*** 

c*** 

c 

c 

c*** 

c*** 

c 

c*** 

c 


c 

Q*** 

c 


COMPUTE  DEL(I)  AND  DELY(I): 

IF(YCEIL(I).GT.YREF(I))THEN 

DEL(l)sMIN((YCEIL(l)-YREF(l))/2.D0,D/2.D0) 

DELY(I) =DABS(YLAY(I)-YVENT) 

ELSE 

DEL(I)=D/2.D0 

DELY(l)s:0.D0 

ENDIF 

COMPUTE  EFFECTIVE  NEAR-VENT  PROPERTIES: 

IF(DELY(I).GT.DEL{I))THEN 

IF(I.EQ.1)THEN 

DEN(1)=DENL(1) 

T(1)=TL(1) 

DO  24  K=1,NPROD 

C{K,1)=CONL(K,1) 

CONTINUE 

ELSE 

DEN(2)<=DENU(2) 

T(2)=TU(2) 

DO  26  K=1,NPROD 

C(K,2)=CONU(K,2) 

CONTINUE 

ENDIF 

ELSE 

IF(I.EQ.1)THEN 

DEN(1)sDENU(1)*(1.DO-(PELY(1)/DEL(1))^APOWER)) 

+ DENL(1  )*(PELY(1  )/DEL(1  ))**APOWER) 

T(1)=TU(1)*DENU(1)/DEN(1) 

DO  28  K=1,NPROD 
C(K,1)-(CONU(K.1)*DENU(1)* 

(1.DO-(PELY(1)/DEL(1))**APOWER)) 

•fCONL(K.1)*DENL(1)*(PELY(1)/DEL(1))**APOWER))/DEN(1) 

CONTINUE 

ELSE 

DEN(2)  3:DENL(2)*(1  .DO-((DELY(2)/DEL(2))**APOWER)) 

+DENU(2)*((DELY(2)/DEL(2))**APOWER) 

T(2)  «TL(2)*DENL(2)/DEN(2) 

DO  30  K-1,NPROD 

C(K^)  « (CONL(K^)*DENL(2) 

*(1  .DO-((DELY(2)/DEL(2))**APOWER)) 
•l■CONU(K^)*DENU(2)*((DELY(2)/DEL(2))**APOWER))/DEN(2) 

CONTINUE 

ENDIF 

ENDIF 

CONTINUE 

DELOEN  IS  DENSITY  IMMEDIATELY  ABOVE  THE  VENT  LESS  DENSITY 
DENSITY  IMMEDIATELY  BELOW  THE  VENT 

DELOEN  « DEN(1)  • DEN(2) 

3.  CALCULATE  VST(0,  THE  "STANDARD"  VOLUME  RATE  OF  FLOW 
THROUGH  THE  VENT  INTO  SPACE  I 

CALCULATE  VST(I)  IF  DELP  « 0 

IFPELP.EQ.O.OOO)THEN 
VST(1)  * 0.0D0 
VST(2)  - 0.000 

ENDIF 

IFPELP.EQ.0.D0)  GOTO  32 
CALCULATE  VST(0  FOR  NONZERO  DELP 
IF(DELP.QT.O.OOO)THEN 
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VST(2)  s O.ODO 
RHO  > DEN  (2) 

EPS  * DELP/P(2) 

ENDIF 

IF(DELP.LTO.ODO)THEN 

VST(1)  = O.ODO 
RHO  = DEN(1) 

EPS  = -DELP/P(1) 

ENDIF 

X = 1.0DO-  EPS 

COEF  « 0.60D0  -t-  0.25D0*EPS 

XO  «:  1.0D0 

EPSCUT  = EPSP*MAX(XO,DPREF(1),DPREF(2)) 

EPSCUT  B DSQRT(EPSCUT) 

SRDELP  B DSQRT(DABS(DELP)) 

FNOISE  B 1.0D0 

IF((SRDELP/EPSCUT).LE.130.D0)THEN 

FNOISE  B 1.0D0  • DEXP(-SROELP/EPSCUT) 

C 

NOTE:  IF  SINGLE  PRECISION  THEN  USED  65.  INSTEAD  OF  130. 

C 

ENDIF 

IF(EPS.UE.0.1D-5)THEN 

W B 1.0D0  • 0.7SD0*EPS/QAM 

ELSE 

IF(EPS.LT.QAMCUT)THEN 

GO  B X**(1.0D0/GAM) 

FF  B DSQRT((2.0D0*GAM/(GAM*1.0D0))*GG*GG* 

$ (1.0D0-X/GG)) 

ELSE 

FF  B GAMMAX 

ENDIF 

W B FF/DSQRT(EPS-I-EPS) 

ENDIF 

V « FNOISE*COEF*W*DSQRT(2.0DO/RHO)*AVENT*SRDELP 
IF(DELP.GT.0.0D0)VST(1)  > V 
IF(DELP.LT0.0D0)VST(2)  b V 
32  CONTINUE 
C 

C***  4.  WHEN  CROSS-VENT  DENSITY  CONFIGURATION  IS  UNSTABLE,  I.E., 

C***  DELDEN  > 0.  THEN  CALCULATE  THE  VENT  FLOW  ACCORDING  TO: 

C***  COOPER,  LY,  COMBINED  PRESSURE-  AND  BUOYANCY-DRIVEN  FLOW 

C***  THROUGH  A HORIZONTAL  VENT,  TO  APPEAR  AS  NISTIR,  NATIONAL 

C***  INSTITUTE  OF  STANDARDS  AND  TECHNOLOGY  GAITHERSBURG  MD 

C 

C***  FOR  STABLE  CONFIGURATION,  GO  TO  S.  AND  VOLUME  FLOW  RATES  WITH 
C***  STANDARD  MODEL 

C 

IF(DELOEN.LE.O.ODO)  GOTO  35 
C 

C***  FOR  UNSTABLE  CONFIGURATION,  NOW  CALCULATE  THE  COMBINED  PRESSURE- 

C***  AND  BUOYANCY-DRIVEN  VOLUME  VENT  FLOW  RATES  FROM  THE  HIGH-TO-LOW 

AND  FROM  THE  LOW-TO-HIGH  SIDES  OF  THE  VENT,  VBHIGH  AND  VBLOW, 

C***  RESPECTIVELY  FROM  THESE  THEN  CALCULATE  THE  COMBINED  PRESSURE- 

C***  AND  BUOYANCY-DRIVEN  VOLUME  FLOW  RATES,  VB(0,  THROUGH  THE  VENT 

C***  INTO  SPACE  I. 

C 

TBARB(T(1)-i>T(2))/2.D0 

DENBARbPEN(1)<I-DEN(2))/2.D0 

XMEW-(0.04128O-7)*OENBAR*(TBAR**2.5D0)/(TBAR-t-110.4D0) 

EPSDENbDELOEN/DENBAR 

ORb2.DO*0*(D**3)*EPSDEN/((XMEW/DENBAR)**2.DO) 

IF(DELP.OT.OOO)  EPSDEN—EPSDEN 

VHIGHFLb0.17S4O0*AVENT*DSQRT(2.D0*O*D*DABS(EPSDEN))* 

1 DEXP(0.553600*EPSDEN) 

DELPFLb0.2427D0*(4.D0*O*DABS(EPSDEN*DENBAR)*D)*(1.D0-I-EPSDEN/2.D0)* 

1 DEXP(1.107200*EPSDEN) 

DPDDPFL«DABS(DELP)/DELPFL 
SIGMAI  bFNOISE*COEF*W/0.1780D0 
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IPPPDDPFLQE.I  .DO)THEN 
VBLOW=O.DO 

VBHIOH  * VHIGHFL*(1  .DO-(SIGM  A2**2.D0)  + 

1 DSQRT(SIGMA2*M.D0+(SIGMA1**2.D0)*  PPDDPFL-I.DO))) 

ELSE 

VEXMAX=0.055D0*(4.D0/PI)*AVENT*DSQRT(G*D*DABS(EPSDEN)) 

XM  « (SIGMA1  /SIGM  A2)**2.D0-1  .DO 

VBLOW=VEXMAX*(((1  .D0+XM3/2.D0)*((1  .DO-DPDDPFL)**2)- 
1 (2.D0+XM3/2.D0)*(1.D0-DPDDPFL))**2.D0) 

IFPPDDPFLEQ.0.D0)THEN 
VBHIGH-VBLOW 

ELSE 

VBHIGH>sVBLOW-t- 

1 (XM-DSQRTO  .DO + (XM**2.D0-1  .D0)*(1  .DO-DPDDPFL)))* 

1 

VHIGHFU(XM-1.DO) 

ENDIF 

ENDIF 

IF(DELP.GT.O.DO)THEN 

VB(1)sVBHiGH 

VB(2)=VBLOW 

ELSE 

VB(1)«VBLOW 

VB(2)sVBHIGH 

ENDIF 

C 

C***  5.  CALCULATE  WENT(I),  THE  VOLUME  RATE  OF  FLOW  THROUGH 

C***  THE  VENT  INTO  SPACE  I 

C 

35  CONTINUE 

IF(DELDEN.L£.O.ODO)THEN 

VVENT(I)  = VST(1) 

VVENT(2)  - VST(2) 

ELSE 

IF((O.DO.LT.OR)AND.(OR.LT.2.07))THEN 
VVENT(I)  «VST(1) + (OR/2.D7)*(VB(1)-VST(1)) 

VVENT(2)  ■ VST(2) + (OR/2.D7)*(VB(2)-VST(2)) 

ELSE 

VVENT(1)«VB(1) 

VVENT(2)cVB(2) 

ENDIF 

ENDIF 

C***  6.  CALCULATE  THE  VENT  FLOW  PROPERTIES 

C 

DENVNT(I)  « DEN(2)*P(1)/P(2) 

DENVNT(2)  > DEN(1)*P(2)/P(1) 

TVENT(I)  ■ T(2) 

TVENT(2)  - T(1) 

DO  50  K - 1,  NPROD 

CVENT(K,1)  - C(K^) 

CVENT(K^)  - C(K.1) 

50  CONTINUE 

C 

C***  7.  CALCULATE  THE  VENT  MASS  FLOW  RATES 

C 

DOM  I « 

XMVENTO)  > DENVNT(I)*VVENT(I) 

M CONTINUE 

C 

C^  8.  CALCULATE  THE  REST  OF  THE  VENT  FLOW  RATES 

C 

DO  70  1-1^ 

QVENTO)  ■ XMVENT(I)*CP*7VENT(I) 

DO  65  K - 1,  NPROD 

PVENT(K,0  ■ XMVENT(l)*CVENT(K,0 
CONTINUE 
CONTINUE 

C***  9.  CALCULATE  THE  RATE  AT  WHICH  THE  VENT  FLOWS  ADD  MASS. 
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c*** 

Q*** 

c*** 

c 


72 


75 


80 


85 


00 


# 


05 


100 


ENTHALPY.  AND  PRODUCTS  TO  THE  LAYERS  OF  THE  SPACES  FROM 
WHICH  THEY  ARE  EXTRACTED.  FIRST  TREAT  LAYERS  OF  SPACE  1 
AND  THEN  SPACE  2. 

DO  150  I s 

FRACM(I)=1.D0 
FRACO(l)«1.D0 
DO  72  KbI.NPROD 

FRACP(K.I)=1.D0 

CONTINUE 

IF(DELY(I).QT.DEL(I))THEN 

IF(I.EQ.1)THEN 

XML(1)«-XMVENT(2) 

XMU(1)«0.0D0 
QL(1)  - ^VENT(2) 

QU(1)  s 0.0D0 
DO  75  Kel,NPROD 

PL(K,1)  = -PVENT(K.2) 

PU(K,1)  > O.ODO 
CONTINUE 

IF(YCEIL(1).EQ.YREF(1))THEN 

XMU(1)sXML(1) 

aU(1)=QL(1) 

DO  80  K=1,NPROD 
PU(K,1)=PL(K,1) 

CONTINUE 

ENDIF 

ELSE 

XMU(2)s-XMVENT(1) 

XML(2)«0.D0 

QU(2)«-QVENT(1) 

QL(2)-0.D0 
DO  85  Ke1.NPROD 

PU(K.2)— PVENT(K.1) 

PL(K4I)«0.D0 

CONTINUE 

IF(YCEIL(2).EQ.YREF(2))THEN 

XML(2)>>XMU(2) 

QL(2)bQU(2) 

DO  00  Ks1,NPROD 

PL(K.2)-PU(K.2) 

CONTINUE 

ENDIF 

ENDIF 

ELSE 

IF0.EQ.1)THEN 

XMU(1)»DENU(1)*(1.D0>((DELY(1)/DEL(1))**APOWER)) 

*VVENT(2) 

XML(1)— DENL(1)*((PELY(1)/DEL(1))**AP0WER))*VVENT(2) 
IF(XMVENT(2).QT.0.D0)FRACM(2)«XML(1)/XMVENT(2) 

QUO)  >XMU(1)*CP*TU(1) 

QL(1)  a>XML(1)*CP*TL0) 

IF(QVENT(2).GT.0.D0)FRACQ(2)>iQL0)/QVENT(2) 

DO  05  K-1,NPROD 

PU(K.1) -XMU0)*CONU(K.1) 
PL<K.1)«XML0)*CONL(K,1) 

IF(PVENT(K.2).QT.0.D0)FRACP(K.2)  «PL(K,1)/PVENT(K.2) 
CONTINUE 

IF(YCEIL0)-EQ.YREF(1))THEN 

XML0)«XMU0) 

QL0)>QU0) 

DO  100  K-1.NPROD 

PL(K.1)-PU(K.1) 

CONTINUE 

ENDIF 

ELSE 

XML(2)«- DENL(2)*0-DO-(PELY(2)/DEL(2))**APOWER)) 
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# 

*WENT(1) 

XiyiU(2)=-DENU(2)*(((DELY(2)/DEL(2))**APOWER))*VVENT(1) 

IF(XM  VENT(1  ).QT.0.D0)  FRACM  (1 ) =XMU  (2)/XMVENT(1 ) 

QL(2) =XML(2)*CP*TL(2) 

QU(2)  =XMU(2)*CP*TU(2) 
IF(OVENT(1).OT.O.DO)FRACQ(1)=QU(2)/QVENT(1) 

DO  105  K=1,NPROD 

105 

PL(K,2)  =XML(2)*CONL(K,2) 

PU  (K^)  sXM  U (2)  *CONU  (K,2) 

IF(PVENT(K,1).OT.O.DO)FRACP(K,1)=PU(K^)/PVENT(K.1) 

CONTINUE 

110 

IF(YCEIL(2).EQ.YREF(2))THEN 

XMU(2)sXML(2) 

QU(2)sQL(2) 

DO  110  K=1.NPROD 

PU(K^)=PL(K^) 

CONTINUE 

ENDIF 

ENDIF 

ENDIF 

150  CONTINUE 

200  CONTINUE 

RETURN 

END 
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